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ABSTRACT 

The  niimerical  weather  analysis  procedures  of  the  National 
Meteorological  Center  are  reviewed,  and  problems  in  adapting  them  to  a 
parallel  computer  are  considered.   Algorithms  for  storage  allocation, 
processing  teletype  data,  dynamically  arranging  data,  and  producing  the 
analyzed  field  are  presented. 


Ill 
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1.   INTRODUCTION 

Niomerical  weather  analysis  was  developed  to  provide  improved  input 
to  numerical  forecasting  programs.  Meteorological  data  is  obtained  by  observa- 
tion from  stations  that  are  located  neither  in  a  regular  pattern  nor  uniform 
distribution.   This  data,  which  is  transmitted  via  teletype  in  standard  code 
formats,  must  be  collected,  decoded,  and  reduced  to  the  grid  used  by  the  fore- 
casting routine.   The  National  Meteorological  Center  has  developed  effective 
procedures  that  produce  operational  analysis  charts,  as  well  as  input  to  fore- 
casting programs. 

The  ILLIAC  IV  computer  will  be  an  array  of  arithmetic  units  operating 
together  \ander  central  control.  In  a  one  quadrant  array,  sixty-four  processing 
elements  can  each  be  doing  two  32-bit  arithmetic  operations  simultaneously. 

Implementation  of  a  numericaJ.  weather  analysis  program  on  ILLIAC  IV 
requires  consideration  of  operation  on  teletype  reports  in  some  parallel  manner, 
storage  of  data  so  that  it  is  readily  accessible  to  appropriate  processing  ele- 
ments, and  maximization  of  simultaneous  operations  across  the  grid. 
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2.      NUMERICAL  WEATHER  ANALYSIS  i 

ij 

2 . 1  Relation  to  Numerical  Forecasting  |i 

From  the  time  von  Neuman  and  Charney  investigated  automated  forecast    i 
procedures,  meteorologists  realized  that  input  data  was  a  major  problem.   Indeed,  j 
part  of  Richardson's  difficulties  in  his  early  attempt  at  manuaJ.  numerical  inte 
gration  to  produce  prognostic  information  was  the  lack  of  upper  air  data   [IJ. 
It  is  the  atmospheric  layers  above  the  surface  to  which  dynamic  and  thermody- 
namic equations  have  been  most  successfully  applied;  upper  air  observing  sta- 
tions are  considerably  fewer  in  number  than  surface  ones.  j 

Manual  analysis  is  highly  subjective  and  time  consuming.   The  same 
objections  apply  to  manual  extraction  of  data  from  charts  for  input  to  machine 
analysis.   In  sharp  contrast,  the  operation  of  a  numerical  forecast  program  is 
quite  fast  and  as  objective  as  the  atmospheric  model  it  simulates. 

2.2  Dynamics  and  Thermodynamics  for  Upper  Air  Analysis 

The  equations  of  motion  applicable  to  a  parcel  of  air  on  the  rotating 
earth  can  be  simplified  by  assuming  that  the  only  forces  present  are  those  due 
to  pressure  gradient,  gravity,  and  the  earth's  rotation,  and  that  motion  reaches 
an  equilibriiun  state  [2].   The  simplified  equations  are 

^    1  ^P      ^     1  ^p  1  ^p  .  ^  ^ 

fv=-r^      fu=--«N^     g=-— =r^+2n  cos  cp 
p  ox  p  oy     *=     p  dz  ' 

v/h-^re  X  is  positive  in  the  eastward  direction,  y  northward,  and  z  upward;  u  and 


Additionally,  Richardson's  work  failed  because  the  compensation  of 
horizontal  forces  was  not  considered  and  because  of  computational  instability 
in  his  model. 


V  are  the  x  and  y  components  of  velocity,  respectively;  p  is  pressure,  p  is 
density,  g  is  acceleration  of  gravity,  and  f  =  2  fi  sin  op  is  the  Coriolis  para- 
meter at  latitude  q?  for  the  earth's  angular  velocity  fi.   The  third  equation  is 

■X- 

the  equation  for  hydrostatic  equilibrium  after  the  small  vertical  Coriolis  term 
is  neglected.   The  horizontal  equations  are  termed  the  geostrophic  wind  equa- 
tions and  express  the  horizontal  balance  of  Coriolis  amd  pressure  gradient 
forces. 

When  the  vertical  coordinate  is  transformed  to  pressure,  the  horizon- 
tal equations  sire 


V  =  f  If  1     and     u  =  -  I 


In  analysis  work  these  equations  have  the  advantage  that  density  is  not  explicitly 
involved.   The  constant  pressure  surfaces  are  sufficiently  quasi -horizontal 
that  they  are  easily  usable  in  place  of  constant  height  surfaces. 
Ihe  hydrostatic  equation,  ignoring  horizontal  forces,  is 

dp  =  -0  g  dz. 

Substituting  the  ideal  gas  equation  of  state  and  integrating  from  z,  to  z^  yields 
[3]  the  hypsometric  equation 

z^  -  z,  =  —  (In  p^  -  In  p_), 
2    1   mg      1      j'2'-' 

where  m  is  the  average  molecular  weight  of  dry  air,  R  is  the  gas  constant,  and 


Ifydrostatic  equilibrium  is  one  of  the  most  basic  assumptions  made  in 
theoretical  meteorology.  It  quantifies  the  fact  that  the  vertical  scale  of  mo- 
tion is  orders  of  magnitude  smaller  than  the  horizontal  one. 


T  is  the  mean  virtual  temperature  of  the  layer.  Virtual  temperature,  T  ,  is  the 
temperature  of  dry  air  of  the  same  density  and  pressure  as  the  given  moist 
sample.   It  rarely  exceeds  the  actual  temperature  by  more  than  three  degrees 
Celcius. 

Thus  the  thickness  of  a  layer  of  air  between  two  constant  pressure  sur- 
faces is  directly  proportional  to  the  mean  virtual  temperature  of  the  layer. 
From  a  plot  of  temperature  and  moisture  content  versus  logarithm  of  pressure, 
the  virtual  temperature  plot  can  be  easily  developed.   Equal  area  averaging  of 
the  virtual  temperature  between  two  pressures  then  yields  the  thickness  of  the 
layer.   Because  of  the  averaging  nature  of  the  hypsometric  equation,  relatively 
simple  upper  air  sounding  equipment  will  yield  good  measurements  of  the  heights 
of  constant  pressure  surfaces. 

The  hydrostatic,  geostrophic,  and  hypsometric  equations  contain  the 
factor  g,  the  earth's  acceleration  of  gravity.  Meteorologists  usually  combine 
the  acceleration  due  to  gravitational  forces  with  the  centrifugal  force  of  the 
earth's  rotation  into  g.  Thus  g  varies  slightly  with  latitude  and  surfaces  of 
equal  geopotential  (specific  potential  energy  due  to  gravity)  do  not  correspond 
to  those  of  geometric  altitude.  The  above  equations  are  usually  formulated  in 
terms  of  geopotential  rather  than  geometric  height.  The  geopotential  meter  is 
defined  from 


'  /o  ^^' 
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and  gives  geopotential  height  measurements  numerically  similar  to  geometric 
heicht  in  meterr,  fP].   Actually,  geopotential  meters  are  a  measure  of  specific 
energy. 


In  much  of  the  development  of  theoretical  meteorology,  specific  vol- 
tune,  specific  force,  etc.,  are  used  and  explicit  measure  of  the  mass  of  air  and 
air  parcels  is  not  required. 


In  sunmiary,  the  hydrostatic  and  geostrophic  equations  express  the 
equilibrium  that  is  the  general  nature  of  the  atmosphere.   The  geostropic 
equations  yield  a  simple  relationship  between  vrind  velocity  and  the  gradient 
of  pressure,  or  equivalently,  of  geopotential  height.   The  hypsometric  equation 
provides  a  simple  means  of  determining  geopotential  height  of  constant  pressure 
surfaces. 

2.3  The  National  Meteorological  Center  (NMC)  Procedures 

Cressman  [h],  later  at  the  Joint  N-umerical  Weather  Prediction  Unit  (the  or- 
ganizational predecessor  of  NMC),  considered  an  analysis  procedure  for  geopoten- 
tial height  based  on  quadratic  least  squares  interpolation  of  observed  data  at 
grid  points.   The  asymmetry  of  data  distribution,  particularly  over  ocean  areas, 
caused  poor  results  in  such  areas  of  sparse  data.  Furthermore,  extensive  compu- 
tation was  required  to  invert  the  matrices  of  coefficients  during  the  computa- 
tion at  each  grid  point. 

Bergthorsson  and  Doos  [5]  developed  a  procedure  at  the  University  of 
Stockholm  based  on  modification  of  a  first -guess  field.   Their  reasoning  was 
that  interpolation  of  any  nature  will  fail  in  an  area  lacking  observational  data 
and  that  in  such  areas  a  forecast  based  on  old  data  is  better  than  an  inter- 
polated analysis.   In  tropical  areas,  climatological  normals  are  relatively  good 
first  guesses. 

Their  procedure  begins  with  determination  of  a  first-guess  field  as 
a  weighted  average  of  the  climatological  normals  and  the  numerical  forecast 
based  on  the  data  from  the  preceding  observation  time.  Next,  three  types  of 
corrections  to  the  first-guess  field  are  determined  for  each  reporting  station. 
The  first  correction  assumes  that  the  difference  at  the  station  of  observed 
height  (z  )  and  first-guess  height  (z  ,  as  determined  by  interpolation)  is 
equal  to  the  difference  at  s\irrounding  grid  points  of  the  actual  height  (z^) 


and  first-guess  height  (Z„) 


^1  =  ^f  ^  ("O  -  "fs^- 


By  assimiing  the  observed  wind  is  geostrophic  and  representative  of 
the  interval  between  the  station  and  the  pointy  a  second  corrected  value  (z  ) 
is  obtained 


\-\*  tejg  "^  '  l57Jo  ^ 


where 


d  =      \/(Ak)2  +    {^f 


is  the  distance  in  grid  units  from  station  to  grid  point  and  (^j  and  [■t~] 
are  the  x  and  y  components  of  height  gradient  as  computed  from  the  observed 
wind. 

The  third  corrected  value  is  similar  to  the  second  but  uses  the  gra- 
dient at  the  grid  point.   The  relative  weighting  (w^  \t  ,   w  )  of  the  three 
corrected  values  was  determined  by  statistical  regression.   For  stations  report- 
ing only  height,  the  correction  z,  is  computed  but  a  different  weighting  func- 
tion is  used. 

Because  of  the  correlation  of  data  from  stations  in  areas  of  high 
station  density,  the  composite  corrected  value  is  further  weighted  by  the  in- 
verse of  station  density.   The  computed  value  at  each  grid  point  is  then  the 
weighted  average  of  all  composite  corrections  for  stations  and  of  the  first- 
guess  value.   One  scan  of  the  grid  was  used  by  Bergthorsson  and  Doos  in  most  of 
their  experimental  work. 


Cressman  [6]  modified  the  Bergthorsson  emd  Doos  procedure  by  adopting 
a  simplified  weighting  function  and  making  multiple  scans  of  the  grid  (see 
Figure  1  for  NMC  grid).   This  modification  also  drops  the  third  type  of  cor- 
rection.  The  weighting  function,  the  same  for  both  types  of  corrections,  is 


2    2 

-  R  -  d 
w  - 

R  +  d 


where  R,  the  radius  of  influence,  is  the  maximum  d  for  which  a  station  is  used 
to  correct  a  grid  point.   This  type  of  weighting  function  effects  a  smoothing 
of  all  data  within  distance  R.   The  value  of  R  is  decreased  on  each  of  four 
consecutive  scans  in  order  that  the  phenomenon  of  decreasing  scale  can  be  in- 
cluded in  the  resulting  final  anal.ysis  while  errors  in  the  first-guess  field 
are  quickly  smoothed  out. 

In  order  that  height  data  may  be  given  greater  weight,  wind  data 
(correction  2)  is  not  included  on  the  first  scan  of  a  height  analysis.   Sta- 
tions reporting  wind  only  are  utilized  by  forming  a  type  2  correction  (z')  with 

z^  replaced  by  z _  . 
0  IS 

Since  196^,  the  value  R  (the  radius  of  influence)  at  a  grid  point  is 
determined  by  data  density  [7].  This  variation  of  R  is  a  discrete  function  of 
the  number  of  height  reports  (see  Table  l). 


TABLE  1 

EXAMPLE  OF  RADIUS  OF  INFLUENCE 
AS  A  FUNCTION  OF  DATA  DENSITY 

R  in  grid  units  Number  of  height  reports 

2.2  5  or  more 

3.5  7  or  more 

6.9  1  or  more 
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The  calculated  value  at  the  grid  point  is 


i  m        n 

AZw.c.  +B  Z-W.C.+  Sw,  c, 
i=l  "  "     ,1=1  J  J   k^l  ^  ^ 

^   ^f  ^  Ai  +  Bm  +  n 


"i  =  ^\  -  "f^i  ^  "O  -  "fs 


^d  =  ^^2  -  ^P,  ^  ^fs  ^  f4"^rlr),^r  ^f 


The  subscripts  i,  j,  and  k  indicate  values  for  one  of  the  i  +  m  +  n  stations 
that  can  influence  the  report.   The  coefficient  A  can  be  regarded  as  the  weight 
given  to  gradients  of  the  first-guess  field  in  relation  to  the  observed  winds, 
weighted  as  1  [6].   A  equals  t^.      B  similarly  is  the  weight  given  to  the  inter- 
polated first-guess  value  as  a  basis  for  extrapolating  the  gradient  [7]. 

B  =  0,  7T,  -r,   and  —  on  scans  1,  2,   3^  and  h,   respectively. 

o  z     o  z 
The  value  of  the  components  of  height  gradient  v-  and  ^^  are  computed 

relative  to  the  grid  coordinates.   If  the  geostrophic  gradients  are  determined 

by 


\^/     g  ^        \5y/     g  ^ 


then  an  approximation  of  non-geostrophic  affects  can  be  included  by  computing 


if  =  ¥1  ^  fi  (-(I)  °i.  ^  (ll>) 


g     I  o        g  g 


Again  f  and  g  are  the  Coriolis  parameter  and  acceleration  of  gravity,  respec- 
tively; m  and  6  relate  to  the  map  on  which  the  grid  is  based.  K  is  a  constant 
less  than  one  which  is  assigned  to  a  constant  pressure  surface  for  a  season. 

D'  ,  D'  ,  and  D'   are  the  second  derivatives  of  the  height  field.   On  the  square 
xx'  yy'      xy  oh 

grid  used  by  NMC,  the  neaj*est  eight  points  to  a  given  grid  point  are  required 
to  compute  these  derivatives. 

Before  each  scan  of  the  grid,  each  station  report  is  checked  for  its 
consistency  with  its  neighbors.   The  mean  difference  of  each  station  within 
two  grid  Tinits  from  the  first-guess  field  is  compared  with  the  difference  for 
the  station  under  consideration.   The  station's  report  is  marked  for  acceptance 
or  rejection  on  the  basis  of  this  comparison  and  data  density. 

The  source  first-guess  fields  used  by  NMC  vary.   The  1000  mb,  500  mb, 
and  200  mb  first-guess  fields  for  height  are  obtained  by  adjusting  in  the  tropic- 
aJ.  regions  the  preceding  forecast  by  climatological  height  gradients.   For  some 
levels,  vertical  extrapolation  or  interpolation  from  levels  previously  analyzed 
provides  the  first-guess  field. 

Except  for  the  500  mb  level,  temperature  is  analyzed  at  the  same  time 
as  height.   The  temperature  analysis  is  like  that  above  but  omits  gradients. 
The  u  and  v  components  of  wind  are  analyzed  by  taking  the  geostrophic  winds  from 
the  corresponding  height  analysis  as  the  first-guess.   Two  scans  are  used  to 
adjust  the  field  to  reported  data.   Other  modifications  to  the  basic  procedure 
are  used  to  accomplish  objective  tropopause  analysis  and  wind  fields  in  the 
tropics  [8].   The  latter  uses  a  grid  on  a  Mercator  map  projection. 
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2.k     The  MC  Grid 

The  grid  used  "by  NMC  (Figure  l)  for  both  analysis  and  forecast  pro- 
cedures is  obtained  by  drawing  one-half  inch  squares  on  a  1:30^000^000  scale 
polar  stereographic  projection  of  the  northern  hemisphere.   This  projection  is 
true  at  60°  N  latitude.   The  resulting  grid  unit  6  is  38I  kilometers  or  206 
nautical  miles.   The  map  factor,  the  ratio  of  projection  image  distance  to 
distance  on  the  earth,  is 


1  +  sin  60' 
m  = 


1  +  sin  -"p 


The  grid  is  oriented  so  that  the  meridians  10°  E,  80°  W,  170°  ¥, 
100°  E  are  along  grid  lines.   The  grid  is  numbered  so  that  the  i  (x  coordinate) 
and  j  (y  coordinate)  zero  lines  intersect  in  the  Indian  Ocean  and  the  North  Pole 
is  i  =  33^  j  =  26.   The  following  formulas  may  be  used  to  convert  the  latitude 
and  longitude  of  any  location  to  grid  coordinates 


R   cos  I'D 

e 
r  = 


1  +  sin  i-p 

i  =  ip  +  r  sin(^  +  lOO) 

j  =  jp  +  r  sin(^  +  100) 

R  =  21.20U  grid  units 
e         ^ 

where  ip  and  jp  are  the  grid  coordinates  of  the  pole  and  ^  is  west  longitude. 

The  MC  analysis  procedures  use  the  1977  points  enclosed  by  the  octo- 
gon  depicted  in  Figure  1. 


R  is  the  radius  of  the  earth  on  the  map  projection  image  plane-- 


R  =  6371.2  km  (1  +  cos  30°)/38l  km  per  grid  unit. 
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Figure  1.   The  NMC  Analysis  Grid. 
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3.      THE  ILLIAC   IV  COMPUTER 

3.1  Organization 

The  basic  concept  of  ILLIAC  IV  is  that  the  rate  of  calculation  on 
a  digital  computer  can  be  increased  by  executing  more  than  one  arithmetic 
operation  simultaneously.   Furthermore^  the  parallel  execution  of  arithmetic 
operations  does  not  require  duplication  of  control  functions.   To  this  end, 
ILLIAC  IV  has  6U  arithmetic  units,  or  processing  elements  (PEs),  under  the 
control  of  one  control  unit  (CU).   The  combination  of  one  CU  and  6h   PEs  is 
termed  a  quadrant  and  one,  two,  or  four  quadrants  can  be  united  to  operate  as 
an  array  [9]. 

3.2  Control  Unit 

The  control  unit  (Figure  2)  fetches  and  decodes  all  instructions. 
Only  address  arithmetic  (2U-bit,  fixed  point)  can  be  done  in  the  CU.   Instruc- 
tions that  only  involve  the  CU  include  branching  instructions  that  allow  the 
entire  quadrant  to  execute  from  another  part  of  the  instruction  stream. 

The  advanced  station  (ADVAST)  in  the  CU  does  the  initial  decoding  of 
instructions.   It  contains  four  accumulator  registers  (ACARs)  and  6k   words  of 
data  buffer  (ADB).   Transfer  instructions  are  entirely  executed  by  ADVAST. 

The  final  station  (FINST)  completes  the  decoding  of  all  instructions 
requiring  PE  action.   This  includes  data  fetches  to  ADVAST.   A  queue  of  up  to 
eight  instructions  is  maintained  between  ADVAST  and  FINST. 

The  instruction  look-ahead  unit  (ILA)  maintains  a  buffer  of  up  to  128 
instructions  and  fetches  all  instructions  from  memory. 

The  test  and  maintenance  unit  (IMU)  is  the  only  path  to  or  from  out- 
side control.   Loading  IMU  registers  can  affect  initialization  of  the  control 
unit. 
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Figiire  2.      Control  Unit. 


Ik 


3.3  Processing  Elements 

I 
Processing  elements  (PEs^  Figure  3)  execute  arithmetic,  boolean,        I 

shift,  and  data  transfer  instructions  in  lock-step  \inder  sequencing  from  FINST. 

Arithmetic  instructions  can  be  either  6U-  or  32-bit  floating  point  or  k8-   or 

I 

24-bit  fixed  point  (see  Figure  h) .      There  are  also  a  limited  number  of  instruc- 
tions for  operation  on  eight-bit  bytes. 

Each  processing  element  is  assigned  20U8  words  of  storage  (PEM).   The 
term  processing  unit  (PU)  denotes  the  PE  plus  its  storage.  J 

The  addressable  PE  registers  are  the  A  (accumulator)  register,  the  R 
(routing)  register,  the  D  (mode)  register,  the  X  (index)  register,  the  B  regis- 
ter, and  the  S  register.   Register  C  (the  carry  register  in  the  arithmetic  logic)  ' 
is  addressable  in  connection  with  byte  operations  by  only  one  instruction.        | 
Register  X  is  I6  bits  in  length;  register  D  is  eight  bits  (see  3«5)«   All  others 
are  6U  bits  long. 

The  S  register  can  be  used  to  store  operands.   Its  I6  low  order  bits 
can  function  the  same  as  register  X  for  index  operations. 

The  6U-bit  wide  common  data  bus  (CDB)  allows  the  CU  to  broadcast  oper- 
ands to  all  PEs.  H 

Register  X  and  register  S  indexing  may  be  done  on  shift  operations  as    | 
well  as  on  memory  fetches  and  stores. 

3.4  Memory  Addressing 

The  control  ■unit  addresses  memory  such  that  the  low  order  six  bits 
indicate  a  particular  PE.   The  preceding  I6  bits  address  a  row  of  64  words  in 
PE  memory  (five  of  these  allow  for  expansion  from  the  present  2  =  20h8   words 
per  PE).   Indexing  by  an  AGAR  is  by  word;  by  register  X  or  S,  by  row  (low  order 
bits  discarded). 
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Figure  3«   Processing  Element. 
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Figure  ^4.   Data  Word  Formats. 


3.5  Mode  Control 

input  to  storage  and  to  three  of  the  programmer  addressable  PE 
registers  (A,  S,  and  X)  can  be  inhibited  by  mode  setting  operations.   Thus 
during  instruction  execution,  PEs  can  be  "turned  off",  and  different  PEs  can 
in  effect  execute  different  instruction  streams  at  different  times.   Note,  how- 
ever, that  even  if  all  PEs  are  disabled,  FINST  sequences  through  the  instruc- 
tions. The   first  two  bits  of  register  D,  the  E  and  El  bits,  are  the  enable 
bits. 

A  test  on  register  A  or  an  index  test  on  register  X  (or  the  low  l6 
bits  of  register  S)  results  in  setting  the  I  or  J  mode  bits  (in  32-bit  opera- 
tions, the  G  or  H  bits  are  also  set  by  the  results  of  the  tests  on  the  inner 
word  of  register  A). 

The  F  and  Fl  bits  reflect  faults  (exponent  overflow,  for  example). 
The  setting  of  an  F  bit  results  in  an  interrupt  in  the  CU. 

The  SET  instruction  is  available  to  set  any  mode  bit  from  the  results 
of  a  logicalAND  or  of  an  inclusive  OR  of  any  mode  bit  and  either  the  E  or  the  El  bit.  For 
example,  SET  I  J  AND  -El  will  set  bit  I  to  the  result  of  J  and  not  El.  Mode 
bits  may  also  be  set  by  broadcast  from  the  CU.   In  this  situation,  each  bit  of 
an  ACAR  goes  to  the  designated  bit  in  the  PE  of  corresponding  number  (i.e.,  bit 
0  to  PE  0,  bit  1  to  PE  1,  etc.). 

3.6  Routing 

Because  a  given  memory  location  is  directly  accessible  to  only  one 
PE,  routing  operations  are  implemented  to  allow  transmissions  of  data  from  the 
designated  register  to  register  R  of  another  PE.  A  route  of  length  one  or 
eight  takes  one  clock  cycle  [10].  Routes  of  other  distances  are  accomplished 
by  a  sequence  of  routes  of  distance  eight  and  one. 
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3.7  Input/Output  and  Array  Control 

The  only  means  of  input/output  for  ILLIAC  IV  is  through  a  disk  sys- 
tem with  a  transfer  rate  of  0.5  billion  bits  per  second  and  maximum  access 
time  of  Uo  milliseconds  [9].   A  Burroughs  65OO  computer  channels  data  to  and 
from  the  disk  system. 

The  B650O  also  has  executive  control  of  the  array.   Through  the  IMU 
it  can  initiate  a  program,  accept  l/O  requests  from  the  operating  program,  re- 
turn results  of  l/O  operations,  and  halt  the  array. 
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h.      STORAGE  ALLOCATION 

U.l  Data  Organization  and  Arittimetic  Efficiency 

In  order  that  the  analysis  of  a  field  on  the  NMC  grid  may  be  executed 
most  efficiently,  it  is  necessary  that  the  number  of  operations  executed  in 
each  PE  be  as  nearly  equal  as  possible.   The  analysis  routine  is  based  on  ad- 
justing the  first-guess  field  at  each  grid  point.   Therefore,  allocation  of 
grid  points  equally  is  a  first  step.   Unfortiinately,  reported  data  has  a  dens- 
ity far  from  uniform.   Furthermore,  there  are  reports,  such  as  aircraft  ob- 
servations, that  will  not  occur  in  the  same  location  with  each  analysis.   Thus, 
some  execution  time  procedure  must  be  devised  to  associate  the  reported  data 
and  the  grid  points. 

Use  of  32-bit  fixed  point  (2U-bit  mantissa  only)  arithmetic  will 
maximize  arithmetic  efficiency  while  maintaining  accuracy  consistent  with  ob- 
servational accuracy.   Furthermore,  l6  bits  of  accuracy  in  the  observed  quan- 
tities is  sufficient  if  they  are  translated  against  some  fixed  value.   For  in- 
stance, 500  milibar  heights  expressed  in  whole  meters  with  a  base  of  3500  geo- 
potential  meters  allows  computation  with  11  bit  quantities  ;  3500  to  7595  is 
sufficient  range  for  this  analysis.   Positive  translation  against  some  lower 
bound  is  preferable  to  positive  or  negative  differences  about  a  mean  or  stan- 
dard value  in  that  the  latter  would  require  programmatic  setting  of  the  sign  bit 
after  unpacking  the  data.   The  organization  of  I6  bit  quantities  in  a  block  of 
four  32 -bit  words  is  shown  in  Table  2. 


^'lultiplication  of  two  11 -bit  fixed  point  quantities  results  in  22 
bits  of  product  and  is  within  the  24  bits  allowed. 


20 


TABLE  2 
ORGANIZATION  OF  DATA  FOR  ANALYSIS 


Word  number 

Contents  of  Upper  Half 

0 

i  coordinate 

1 

average  temperature 

2 

observed  temperatiore 

3 

1  wind  component 

Contents  of  Lower  Half 

j  coordinate 
average  height 
observed  height 
j  wind  component 


^.2  Arrangement  of  Grid  Points 

The  assignment  of  grid  points  to  PEs  is  governed  by  two  somewhat  con- 

n     * 
flicting  requirements.   With  1977  points  to  be  assigned  to  12o  PEs  ,   with  four 

quantities  (height,  temperature,  and  two  components  of  height  gradients  as  com- 
puted from  observed  winds)  at  ten  standard  constant  pressure  levels  required 
by  the  analysis  system,  one  report  may  present  forty  quantities  to  be  used  in 
various  analyses.   If  only  200  of  the  better  than  1000  reports  per  analysis 
contain  all  Uo  data  items,  there  would  be  far  more  data  items  than  can  be  con- 
tained in  one  PE.   The  total  report  must  then  be  kept  in  one  location  only. 
Since  a  complete  analysis  at  one  level  of  both  temperature  and  height  can  be 
accomplished  with  only  eight  pieces  of  data  ,  this  data  can  be  extracted  from 
each  report  and  associated  with  all  appropriate  grid  points.  Blocks  of  15  to 
l6  points  must  be  arranged  so  that  on  one  hand  the  total  number  of  stations  to 
be  used  by  the  points  is  evenly  distributed,  and  on  the  other  hand  the  total 


Actually  "half"  PEs,  the  inner/outer  portions  of  the  arithmetic 
lonits  will  be  termed  "PEs"  as  their  functions  during  the  analysis  procedure  are 
essentiaJ-ly  independent.   The  inner  words  will  be  referred  to  with  even  numbers; 
the  outer,  with  odd. 

In  addition  to  the  four  observed  quantities,  station  location  and 
the  first-guess  height  and  temperature  values  (as  interpolated  from  surround- 
ing grid  points)  are  required. 


i 
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number  of  stations  assigned  to  each  PE  is  kept  to  a  reasonable  number.  A 
totally  random  distribution  of  the  points  for  a  PE  across  the  hemisphere  would 
be  disastrous.  With  each  point  possibly  being  influenced  daring  analysis  by 
stations  up  to  6.9  grid  units  away  [7]  such  a  distribution  wo\jld  force  each  PE 
to  store  practically  every  station.   An  arrangement  of  four  by  four  blocks, 
however,  would  result  in  blocks  over  western  and  central  E\arope  with  every  point 
being  influenced  by  many  stations. 

Figure  5  indicates  an  empiracally  determined  arrangement  that  not  only 
baJLances  the  above  two  objectives  but  also  results  in  a  simple  algorithm  to 
interchange  grid  point  values  for  neighboring  points.   The  blocks  are  eight 
points  long  in  the  increasing  i  direction,  two  wide  in  the  increasing  j  direc- 
tion.  The  points  indicated  by  o   in  Figure  5  are  those  of  smallest  i  and  j 
values  for  each  group  of  I6. 

Coverage  of  the  octagon  used  by  NMC  is  not  complete;   ^0  points  at 
boundaries  are  omitted.   I9  of  these  points  are  the  j  =  51  row  and  only  one 
(i  =  49,  J  =  10)  is  as  distant  from  the  octagon  boundary  as  the  diagonal  of  a 
grid  square.   Only  two  rows  of  three  blocks  each  do  not  start  in  columns 
i  =  9,  17,  25,  33^  ^1^  or  k<^,   and  these  six  are  exactly  four  points  "out  of 
phase  ". 

U.3  Neighboring  Points  and  Responsibility  Areas 

With  the  assignment  of  points  for  which  analysis  computations  will  be 
done,  the  problems  of  determining  average  (interpolated)  field  values  and  non- 
geostrophic  height  gradients  for  each  station  can  be  solved.   In  Figure  6,  the 
encircled  points  are  those  to  be  assigned  to  a  PE  for  the  analysis  calculation. 
The  inner  rectangle,  three  by  nine  grid  squares  in  size,  covers  the 
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Figure  5.  Assignment  of  Grid  Points  to  PEs. 
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Figiire  6.  ArrsLngement  of  Points  within  PE. 
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responsibility  area  for  that  PE  .  The  stations  falling  within  this  area  will  be 
the  responsibility  of  the  FE  for  calculation  of  average  values  and  height  gra- 
dients and  for  consistency  check.  An  array  of  length  9I  in  each  FE  provides 
analysis  field  data  for  these  requirements. 

The  consistency  check  of  observed  data  against  first-guess  field 
values  requires  inclusion  of  the  outer  rectangle -points  0  to  13,  25,  26,  38, 
39,  51,  52,  6k,   65,  and  78  "to  90.   The  inclusion  of  these  points  allows  inter- 
polated field  values  to  be  found  for  any  station  within  two  grid  units  of  the 
responsibility  area. 


A  station  falling  on  a  boundary  would  be  assigned  to  the  higher 
number  PE. 
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5.   PROCESSING  TELEIYPE  REPORTS 

5.1  Structure  of  Teletype  Messages 

Figure  7  shows  a  portion  of  a  teletype  message  that  contains  two 
upper  air  reports — those  from  Topeka  (TOP)  and  Dodge  City,  Kansas,  for  1200 
Greenwich  Mean  Time  ((3^),  May  8,  1969*   Appendix  A  gives  a  grammar  for  this 
type  of  report,  the  type  transmitted  by  North  American  stations  [11] ,  inBNFmeta- 
language.  Ihe  syntactic  structure  is  relatively  simple.  Thus  it  is  possible  to  apply 
the  methods  of  He  imerdinger  [12]  to  scan  and  parse  this  language  of  teletype  reports. 

Unfortunately,  however,  this  is  not  the  only  form  for  messages  trans- 
mitted over  the  national  and  international  meteorological  circuits.   Because  of 
the  slow  rate  of  teletype  transmission,  it  would  also  be  undesirable  to  tie  up 
ILLIAC  IV  with  the  task  of  scanning  incoming  teletype  messages.   It  is  therefore 
appropriate  to  use  the  B650O  computer  to  receive  the  messages.  While  sorting 
out  types  of  messages,  the  B6500  can  remove  the  communications  system  portions 
of  the  incoming  stream.   The  first  three  lines  in  Figure  7,  as  well  as  the 
transmitting  station  identifiers  (TOP  and  DDC) ,  are  such  communications  informa- 
tion and  not  part  of  the  standard  meteorological  code.  They  do,  however,  identify 
the  type  and  area  covered  by  the  message  (UJ  is  a  standard  upper  air  report  and 
US  is  self-explanatory).   The  B65OO  would  block  all  reports  of  the  same  type 
together . 

5.2  Scanning  and  Recognition 

Heimerdinger ' s  procedure  checks  for  tokens  of  the  language  by  loca- 
ting delimeters  (in  some  cases  the  delimeters  are  tokens  themselves).   In  the 
language  of  teletype  reports,  delimeters  have  no  other  meaning;  the  tokens  are 
the  five  character  data  groups  plus  the  identifying  TT,  W,  and  QQ  groups.   The 
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NNNN 


ZCZCWNfcLC449C  IAA6/L 
UJUS  KSTL  081410Z 

TOP  TT  5B121  72456  99983  12222  36015  0fcl24  /////  /////  85473 

05203  00531  70040  02738  30018  50565  16930  26531  4072b  29345 

25529  20927  457//  20041  ^0190  471//  24562  15378  533//  25555 
10636  583//  25529 

88231  571//  23072  ' 

77223  24076(5 


VV  5812/  72456  209^3  122-22  1^553  10415  22820  03600  33733 

00100  44711  01505  55696  03149  66645  04245  77572  10313  88499 
16930  99400  29349  11333  40147  22255  543//  33231  571//  4421,8 

^93//  55l98  471//  66170  513//  77132  551//  68120  545//  99100 
583//  31313  25046  549//  20558<p 


OQ  58120  72456  90023  36015  36026  36026  90467  00520  01023 

36C20  9Q^9/  33016  32018  91246  27518  26520  27027  9205/  :27535 
23034  93056  20040  22064  24076  950//  24542© 

DDC  TT  58121  7245  1  99929  07220  01020  00169  /////  /////  85509 
£i3400  03031  '70070  03905  29014  50566  19356  26522  40728  3256! 
26019  3^924  479//  25021  20189  501//  24541  15276  527//  25046 
10636  567//  25520 

80276  515//  24528  "  ■ 

77999® 


VV  5812/  72451  00929  07220  11913  046;^9  22887  02^00  33850 

00400  44824  03000  55742  00500  66701  03700  77680  05561  168669 
04757  99584  19967  11529  15562  22460  21956  1:3400  32561  44352 

39956  55278  515//  66178  497//  771Bf5  567//  31313  V5e43  509// 
24037(L 


GQ  58120  72451  90034  01020  01529  02038  90678  04fe)23  04019 

36010  909//  31012  91246  27516  265  16  28016  9205/  255|?6  26519 
9305/  25522  23544  9426/  2*1037  250^9  950//  250340) 


Figure  7.   Sample  Teletype  Data. 
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reports  are  stored  in  memory  in  increasing  word  addresses  (i.e.,  across  PEs), 
which  is  the  easiest  method  of  organization  for  an  ILLIAC  IV  disk  file  [13] • 
Groups  may  be  split  across  PE  boiindaries.  Each  five  bit  teletype  signal  is 
converted  to  a  corresponding  8-bit  code  and  one  character  is  stored  per  byte. 
With  the  provision  that  delimeters  are  not  tokens  themselves  and  thus  can  be 
disregsirded  after  the  tokens  are  formed,  a  maximum  of  four  tokens  can  be  found 
in  one  8-byte  word. 

The  first  routine  in  Appendix  B  indicates  how  the  scanning  and  recog- 
nition begins.   The  routine  is  coded  in  the  ILLIAC  IV  assembly  language  ASK 
[lU].   At  the  label  NEXTTTYLINE  a  row  of  512  input  characters  is  read  from 
memory.  The   byte  logical  test  operations  are  used  to  build  a  control  string. 
A  reset  bit  seven  in  each  byte  indicates  a  delimeter  present  and  a  reset  bit 
six  indicates  a  delimeter  in  the  preceding  byte.   The  symbol  "®",  which  indi- 
cates the  end  of  an  independent  part  of  meteorological  code,  is  also  searched  for 
and  marked  in  the  control  string. 

At  NEXTBYTE  the  process  of  examining  the  control  string  begins.  The 
control  bytes  are  successively  loaded  into  Register  D  in  order  to  set  the  mode 
pattern.  This  control  of  mode  is  used  to  control  which  PEs  will  be  enabled  to 
-oncatenate  a  new  character  onto  the  old  partial  group  and  which  to  store  com- 
pleted groups  (tokens). 

Any  token  not  completed  at  the  end  of  the  word  (a  non-delimeter  in 
:\e   last  byte)  is  kept  for  possible  concatenation  with  an  incomplete  token  in 
he  next  higher  PE.   After  the  comment  SECTION  JOINER,  each  PE  routes  to  its 
predecessor  the  first  token  it  found  and  the  control  byte.   Only  PEs  with  a 
non-delimeter  last  byte  and  a  routed  token  are  enabled  for  concatenation  of 
"he  token.   A  routed  token  that  is  not  concatenated  (delimeter  last  byte)  or 
hat  is  too  long  (an  error  condition)  is  routed  back  and  stored  in  the  higher 
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The  question  of  efficiency  for  this  routine  is  debatable.   On  one 
hand,  all  PEs  participate  in  the  frequent  tests.   However,  on  a  serial  machine 
transfer  to  another  section  of  code  would  occur  where  disabling  PEs 
occur  on  ILLIAC  IV.   Considering  that  on  a  long  term  average,  PEs  will  be 
executing  store  operations  for  only  three  out  of  every  eight  bytes,  efficiency 
of  20-30'^  might  be  reasonable.   A  time  estimate  for  processing  one  row  in  this 
manner  is  9OO  clock  cycles  or  h^   microseconds  . 

As  a  result  of  the  processing  so  far  described,  five  character  groups 
have  been  recognized  and  stored  left  justified  in  PE  memory.  The  storage  loca- 
tions for  all  PEs  are  a  block  of  four  words  per  PE. 

The  next  step  is  to  place  all  groups  of  one  report  into  a  single  PE's 
memo2:y.   However,  this  process  is  more  efficiently  accomplished  if  most  PEs 
have  found  a  start  of  a  message  group.   Thus  the  recognition  of  lines  continues. 
Approximately  fifty  lines  of  512  characters  can  be  so  processed. 

Although  not  coded,  the  following  algorithm  will  provide  reasonably 
efficient  code.   Relocation  of  groups  is  accomplished  as  follows: 

1)  The  CU  causes  each  PE  to  test  each  group  it  has  located  for  a 
"TT"  group.   The  cummulative  test  i-e suits  are  transmitted  to  the  CU  and  each 
PE  with  a  TT  group  becomes  a  receiver  for  the  remainder  of  that  report. 

2)  The  resulting  pattern  is  used  by  the  CU  to  control  routes  to  the 
PEs  that  contain  TT  groups. 

3)  The  index  of  the  four  word  block  (corresponding  to  one  original 
row  of  teletype  data)  is  routed  to  the  next  higher  PE. 
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All  time  estimates  are  based  on  a  clock  cycle  of  fifty  nanoseconds. 
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k)     The  CU  shifts  the  pattern  of  receiving  PEs  by  one  in  order  to 
locate  the  PEs  that  are  to  transmit  groups;  if  PE  0  is  one  such  PE  (the  route 
was  from  PE  63),  only  it  is  enabled  to  increment  the  fetch  index  by  four  (i.e., 
continue  the  message  on  the  next  row  of  memory). 

5)  The  PEs  with  routed  indices  are  enabled  and  route  the  four  loca- 
tions of  the  indicated  block  to  the  lower  PE.   If  the  receiving  PE  detects  a 
TT  group,  it  has  completely  received  one  report  and  is  disabled  from  further 
stores. 

6)  The  index  is  routed  again  to  the  next  higher  PE  and  the  procedure 
loops  to  step  3  with  a  route  of  2,  then  3^  etc. 

7)  The  CU  periodically  samples  the  modes  of  all  PEs  and  when  all 
receivers  are  disabled,  the  procedure  terminates. 

It  is  estimated  that  this  routing  can  be  done  with  50-70'^  efficiency. 
Based  on  an  average  routing  time  of  four  clock  cycles  and  100  groups  per  report, 
70  to  100  reports  could  be  rearranged  in  this  manner  in  100  microseconds. 

5.3  Decoding 

The  rigid  format  of  one  particular  type  of  report  easily  allows  de- 
coding in  parallel.   After  the  tests  for  a  particular  type  group  (perhaps  for 
the  third  group  of  a  TT  report  being  00  <digit>  <digit>  <digit>)  PEs  that  pass 
the  tests  are  disabled  and  failures  recorded  for  use  by  error  recovery  routines. 
The  successful  PEs  are  then  enabled  and  data  portions  of  the  group  are  converted 
to  binary.   At  this  stage,  data  are  placed  into  32 -bit  words  with  height  and 
temperature  in  one  word  and  the  two  wind  components  in  another. 
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The  determination  of  station  location  is  based  on  directory  compari- 
son.     After  building  in  each  PE  an  index  table  based  on  the  last  digit  of  the 


( 
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station  numbers^  the  CU  broadcasts  the  directory  entries  for  each  of  the  last     t 

digits.   With  128  PEs,  a  typical  set  of  500  to  600  surface  stations  can  be        i 

i 
checked  with  only  one  or  two  comparisons  per  directory  entry.   It  would  be 

indeed  imusual  if  ten  stations  of  the  same  last  digit  were  stored  together  in 

one  PE. 

Once  station  location  is  determined  in  i  and  j  coordinates,  the  wind 
data  for  each  reported  level  can  be  converted  to  the  i  and  j  components  of  geo- 
strophic  gradient.   The  sine  of  latitude,  calculated  to  find  the  geostrophic 
components,  is  saved  to  be  used  in  the  calculations  of  non-geostrophic  gradients 
during  the  analysis  routines. 

As  the  array  of  decoded  data  is  built,  the  data  actually  present  in 
a  report  is  recorded  for  later  use. 

The  basic  recognition  algorithm  should  be  better  than  ^0%  efficient. 
Perhaps  2'^%  of  the  total  time  would  be  spent  waiting  during  memory  fetches  for 
which  no  successful  match  will  occ\ir.  However,  the  directory  lookup  can  be  no 
more  efficient  than  1/128  because  a  station  will  only  be  found  in  one  PE. 

5.^  Hydrostatic  Checking 

The  check  of  the  report  with  the  hydrostatic  formula  can  be  done  in 
parallel  on  128  stations  at  one  time.   However,  as  with  all  errors,  at  best 
invalid  data  can  be  flagged  for  later  processing. 

5.5  Error  Recovery 

The  types  of  errors  possible  are  many  and  the  assumptions  made  con- 
cerning them  will  determine  to  what  degree  parallel  operation  can  be  effected. 
At  best,  only  a  few  PEs  would  be  involved  for  one  particular  type  of  error. 
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6.   RELOCATION  AND  ANALYSIS  OF  DATA 

6.1  Initial  Assignment  of  Stations  to  PEs 

As  the  controlling  routine  of  the  analysis  procedure  sequences  to  a 
new  level,  the  data  stored  in  the  total  report  of  each  station  must  be  ex- 
tracted and  passed  to  all  PEs  with  points  within  6.9  grid  units  of  the  station. 
Hie  location  of  the  station  reports  is  totally  a  function  of  the  original  tele- 
type message  stream.   However,  broadcast  of  station  location  from  the  CU  allows 
128  simultaneous  tests  for  the  location.   The  data  can  then  be  broadcast  with 
only  PEs  passing  the  location  test  enabled  to  store  the  data. 

The  most  efficient  way  to  transmit  data  to  the  CU  is  by  the  indexed 
block  fetch  (BINX)  instruction.   Preliminarily,  each  PE  checks  for  the  first 
station  with  data  for  the  indicated  level.   Stations  with  both  height  and  wind 
information  are  transferred  first,  then  height  only  and  wind  only.   The  source 
of  data  for  the  BINX  instruction  is  a  skewed  array  formed  by  placing  the  loca- 
tion half  word  in  array  location  0  of  the  PE  of  origin,  the  sine  of  latitude  in 
location  1  after  a  route  of  +1,  the  height  and  temperature  in  the  location  2 
after  a  route  of  +2,  and  the  wind  (height  gradient)  components  in  location  3 
after  a  route  of  +3.   The  indices  of  the  PEs  are  set  so  that  one  block  fetch 
•.'ill  transmit  four  complete  sets  of  data  (two  32-bit  words  in  each  location 
fetched). 

As  each  station  is  received  by  a  PE  it  is  further  checked  for  loca- 
*:ion  within  its  area  of  responsibility  or  in  the  neighborhood  area  required  for 
the  consistency  check;  appropriate  index  lists  are  built. 

■•2  Interpolated  Values  and  Consistency  Check 

Before  actual  analysis  can  begin,  the  average  field  values  for  each 
.station  Tirust  be  computed.  The   responsibility  index  list  facilitates  this 
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calculation.   The  sine  of  latitude  is  removed  from  its  location  and  is  re- 

placed  by  the  interpolated  height  and  temperature  values.   The  geostrophic       ; 

wind  components  are  also  extracted  and  the  two  are  saved  in  arrays  ordered 

1 
indentically  to  the  responsibility  list.   The  consistency  check  can  then  use     j 

the  calculated  averages  as  it  compares  for  consistency  of  the  observed  data.      ! 

The  interpolated  field  data  is  then  transmitted  to  all  other  PEs 

using  the  station.   This  transfer  requires  routing  first  the  location  word 

(to  identify  the  station)  and  then  the  interpolated  data.   The  routing  must  be   | 

< 

done  for  the  rows  of  l6-point  blocks  up  to  +k  and  -3  away  and  to  the  adjacent  PEs  of  the   * 

f' 
same  row.  Routing  between  the  inner  half  and  outer  half  of  the  same  PE  is  actually  the 

execution  of  a  SWAPA  instruction^  which  interchanges  the  inner  and  outer  words 
of  the  A  register.   A  sequence  of  routes  and  interchanges  will  cause  the  sta- 
tion report  in  question  to  pass  through  many  PEs  (the  differences  in  lengths 
of  rows  cause  this  situation).   Broadcast  of  a  mode  pattern  controls  the  store 
of  the  resiolts  of  the  distance  test  so  that  the  proper  mode  test  resiilts  can 
be  restored  when  the  interpolated  data  is  routed  and  ready  to  be  stored. 
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6.3  Interchange  of  Field  Values  '  " 

The  first-guess  field  is  initially  loaded  only  into  the  PEs  in  which 
the  calculation  is  to  be  done.   Before  each  scan,  field  values  must  be  trans- 
ferred into  all  PEs  requiring  the  points  for  consistency  check  and  height  gra- 
dient calculation.   The  method  is  essentially  that  of  6.2;  however,  it  is 
programmatic ally  known  from  and  to  which  array  locations  the  data  transfer 
involves.   Thus  no  tests  are  required;  broadcast  of  appropriate  mode  patterns 
will  enable  only  the  desired  stores  after  each  route  and  interchange  operation. 
This  interchange  is  also  simpler  than  that  of  6.2  because  routes  to  adjacent 
rows  are  the  greatest  required. 
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6.U  Analysis  Calculations 

Once  the  transfer  of  data  to  appropriate  PEs  is  complete,  the  cal- 
c\ilation  of  the  analysis  is  straight -forward.  Appendix  B  contains  listings 
of  two  routines  necessary  to  the  actual  analysis. 

The  routine  STATIONSORT  builds  the  array  of  pointers  to  stations  to 
be  used  for  analysis  at  each  point.   Note  that  the  order  of  points  is  not  that 
of  their  occurrence  in  the  array,  but  rather,  an  order  designed  to  improve  com- 
putationetl  efficiency.   That  is,  the  point  within  each  PE's  block  of  l6  for 
which  the  greatest  number  of  stations  are  used  in  the  correction  computation 
is  first  and  the  one  with  the  least  number  of  stations  is  last.   At  the  label 
ALL  in  the  procedure,  one  of  the  foirr  bit  pointers  is  extracted  from  the  two 
word  array  POINTORDER  (see  Figure  8). 

The  first  step  accomplished  with  the  data  for  a  grid  point  is  to 
calcvilate  the  distance  from  grid  point  to  station  (label  AL2).  Thereafter  a 
sequence  of  tests  compares  this  distance  with  each  of  the  three  radii  of  in- 
fluence. For  the  lowest  of  the  radii  of  influence,  RLOW,  and  data  with  both 
height  and  wind  reports,  stations  are  stored  directly  into  the  pointer  array 
STAPERPNT.  For  the  other  combinations  of  radii  and  reported  data,  temporary 
arrays  are  built.  After  all  stations  with  height  reports  are  checked,  wind 
only  reports  are  processed. 

The  array  STAPERPNT  packs  four  8-bit  pointers  to  each  word.   Since 
the  pointers  are  to  four  word  blocks  in  a  PE's  memory,  they  must  be  left  shifted 
before  use.   Also,  in  order  to  use  an  all  zero  pointer  as  a  delimeter  between 
types  of  data,  the  pointers  are  initially  set  one  block  (four  rows  of  PE  memory) 
high. 

The  constant  checking  of  the  same  stations  against  points,  obviously 
leads  to  some  con^jutational  inefficiency.  However,  the  routine  that  assigns 
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stations  to  grid  points  has  a  beneficial  effect  here — all  points  do  not  have 
to  be  checked  against  all  stations.   An  estimate  of  30-50^  efficiency  is  rea- 
sonable. With  the  possibility  some  PEs  have  200  stations,  the  time  for  this 
procedure  approaches  one  fourth  of  a  second. 

The  routine  SCAKFIELD,  however,  is  highly  parallel.   The  procedure 
UNPACK  takes  pointers  to  STAPERPNT  and  locates  the  station  data  into  temporary 
storage.   The  computation  that  follows  is  a  straight -forward  application  of  the 
analysis  formulas.   The  building  of  STAPERPNT  was  such  that  the  computations  of 
the  three  types  of  corrections  are  done  sequentially.   The  final  computed  value 
replaces  the  previous  first -guess  value. 

The  efficiency  of  computation  begins  at  lOO/o  as  every  block  will 
contain  some  data.   However,  it  will  quickly  drop  off.   The  distribution  of 
data  will  probably  force  efficiency  under  '^0%.      The  time  of  one  scan  over  the 
field  is  estimated  to  be  75  milliseconds. 

6.5  Analyst  Interaction 

During  the  execution  of  the  system  presently  employed  at  NMC,  an 
arialyst  monitors  the  machine  analysis.   In  particular,  he  has  the  option  of 
forcing  into  the  analysis  bogus  or  doubtful  data.   Such  procedures  would  be 
wasteful  of  ILLIAC  IV  processing  power.   Considering  the  speed  with  which  one 
scan  can  be  made,  it  is  suggested  that  more  than  one  complete  set  of  analyses 
cam  be  computed,  with  the  analyst  inserting  bogus  data  between  runs  in  order 
to  control  the  results  as  desired. 
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7.   CONCLUSIONS 

The  weather  analysis  problem  is  primarily  one  of  much  data,  assym- 
metrically  distributed,  and  a  few  calculations.   Algorithms  for  storage  of 
unevenly  distributed  data  can  be  developed  for  ILLIAC  IV.   The  handling  of 
input  data  and  calculation  of  the  analysis  can  be  done  with  a  great  deal  of 
parallelism.   The  most  difficult,  and  time-consuming,  steps  are  those  requiring 
broadcast  of  data  from  the  CU.   More  detailed  knowledge  of  the  input  data  may 
lead  to  more  efficient  means  of  directory  look-up  and  of  assignment  of  station 
reports  to  grid  points. 

Considering  that  several  time- con s\;iming  routines  have  not  been  coded, 
an  estimate  for  total  analysis  time  for  one  level  is  five  to  ten  seconds. 
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APPENDIX  A 
SYNTAX  FOR  UPPER  AIR  REPORTS  FOR  LAND  STATIONS 

1.  Report  Structure 

<upper  air  report  for  land  station>  ::=  <upper  level  report  for  land  statlon> 

<line  dellraiter> 

<upper  wind  report  for  land  station>; 
<upper  level  report  for  land  statlon>  ::=  <stajidard  level  report> 

<line  delimeter> 
<significant  level  report>. 

2.  Standard  Level  Report  Structure 

<standard  level  report>  ::=  <identification  section  1>  <delimeter> 

<surface  data>  <delimeter>  <standard  level  data> 
<delimeter>  <tropopause  report>  <delimeter> 
<maximuin  wind  report>  0  ; 
<identification  section  1>  ::=  TT' <delimeter>  <date-time  indicator> 

<wind  levels  indicator>  <delimeter> 
<station  designator> 

<surface  data>  ::=  99<pressure>  <delimeter>  <temperature  group>  <delimeter> 

<wind  group>; 

<standard  level  data>  ::-  <level  and  height  group>  <delimeter> 

<temperature  group>  <delimeter>  <wind  group>; 
<level  and  height  group>  ::=  <level  indicator>  <height>} 
<level  indi.cator>  : :  =  00  1  ^5  |  70  |  50  |  Uo  |  30  |  20  |  15  |  10  ; 
<height^  ::=  <digit>  <digit>  <digit>; 
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<tropopause  report>  ::=  <tropopause  data>  |  <tropopause  report>  <delijneter> 

<tropopause  data>  |  88999 j 
<tropopause  data>  :  :  =  <tropopause  group>  <delijneter>  <temperature  group> 

<delimeter>  <wind  group>; 
<tropopause  group>  ::=  88<pressure>; 
<jnaxim\am  wind  report>  ::=  <top  of  soimding  report>  |  <internal  maxim-um  report>  | 

77999  i 
<top  of  sounding  report>  :  :  =  <top  of  sounding  group>  <delimeter>  <-wind  group>; 
<top  of  sounding  group>  ::=  66<pressure>; 
<intem5Ll  maximum  report>  :  :  =  <internal  maximum  data>  |  <internal  maximum  report> 

<delimeter>  <internal  maximum  data>; 
<internal  maximum  data>  ::=  <internal  maximum  group>  <delimeter>  <wind  group>; 
<internal  maximum  group>  : : =  77<pressure>. 


3.   Significant  Level  Report  Structure 

<significant  level  report>  ::-   <identification  section  2>  <delimeter> 

<significant  temperature  level  report>  <delimeter> 
<250  milibar  report>  <delimeter> 
<regional  code  section>  (D  ; 

<identification  section  2>  ::=  W<delimeter>  <date  time  indicator>/<delimeter> 

<station  designator>; 

<significant  temperature  levels  report>  ::=  <surface  temperature><delimeter> 

<signif leant  temperature  levels  data>; 
<surface  temperature>  ::= 
<significant  temperature  levels  data>  ::=  <signif leant  temperature  level  group> 

<delimeter>  <temperature  group>  | 

<significant  temperature  levels  data> 
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<deliineter> 

<significant  temperature  level  group> 

<dellmeter>  <temperat-ure  group>; 
<signif leant  temperature  level  group>  ::=  <signif leant  level  lndlcator> 

<pressure>; 
<slgnlf leant  level  indicator>  : :=  11  |  22  |  33  |  ^^  |  55  |  66  |  77  |  88  |  99  ; 
<250  mlllbar  report>  ::=  31313<dellmeter>  25<helght>  <dellmeter> 

<temperature  group>  <dellmeter>  <wlnd  group>; 
<regional  code  section>  ::=  51515<dellmeter>  <reglonal  code  groups>; 
<reglonal  code  groups>  : :=  <101  group>  |  <super  remarks>; 
<101  group>  :  :  =  lOKlOl  code>; 
<101  code>  ::=  <diglt>  <digit>; 

<super  remarks>  : : =  SUPER<dellmeter>  <super  layers>; 

<super  layers>  ::=  <dlgit>  <dlglt>  -  <dlglt>  <digit>  |  <super  layers>  <dlgit> 

<dlglt>  -  <dlglt>  <dlgit>. 

h.      Upper  Wind  Report  Structure 

<upper  wind  report  for  land  station>  ::=  <reglonal  standard  level  report>  ®  ; 
<regional  standard  level  report>  ::=  <ldentlflcatlon  section  3>  <dellmeter> 

<reglonal  standard  level  data>; 
<i  dent  ifi  cat  ion  section  3>  ::=  QQ<delljneter>  <date  time  lndlcator> 

<equipment  indicator>  <dellmeter> 
<statlon  deslgnator>; 
<equipinent  indicator>  : :  =  <digit>; 
<regioneLL  standard  level  data>  ::=  <level  groups>  | 

<regional  standard  level  data>  <delimeter> 
<level  groups>j 
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<level  groups>  ::=  <single  level  groups>  |  <double  level  groups>  | 

<triple  level  group s>; 
<single  level  groups>  ::=  9^tens  digit>  <ones  digit>  //  <delimeter>  <wind  group>; 
<double  level  groups>  ::=  9*^tens  digit>  <ones  digit>  <ones  digit>  /  <dellmeter> 

<vrind  group>  <delimeter>  <wind  group>; 
<triple  level  groups>  ::=  9"^tens  digit>  <ones  digit>  <ones  digit>  <ones  digit> 

<deliineter>  <wind  group>  <delimeter>  <wind  group> 

<delinieter>  <wind  group>; 
<tens  digit>  : : =  <digit>; 
<ones  digit>  ::=  <digit>. 

5.  General  Terminal  Symbols 

<delimeter>  : : =  <line  delimeter>  |  <blank>j 

<line  delimeter>  ::=  {combination  of  the  carriage  retiirn  and  line  feed  teletype 

control  symbols  that  cause  typing  to  begin  on  a  new  line} ; 

<blank>  ::=  (teletype  impulse  that  spaces  the  typing  equipment  without  printing 
a  character] . 

6.  General  Non-Terminal  Symbols 

<date-time   indicator>   ::=  <date  and  units  code>  <time>; 

<date  smd  units   code>   ::=  <digit>  <digit> 

<time>   ::=  00    |    12    ; 

<wind  levels  indicator>  ::=  <digit>  ; 

<station  designator>  ::=  <block  number>  <station  number>  ; 

<block  number>  ::=  <digit>  <digit>; 

<station  number>  ::=  <digit>  <digit>  <digit>; 
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<press-ure>  ::=  <dlgit>  <digit>  <digit>j 

<temperature  group>  ::=  <temperature  in  degrees>  <tenths  and  sign  lndicator> 

<de-w  point  depression>  |  /////; 
<temperature  in  degrees>  ::=  <digit>  <digit>  <digit>; 
<tenths  and  sign  indicator>  ::=  <digit>; 
<dew  point  depression>  ::=  <digit>  <digit>; 
<wind  group>  ::=  <wind  direction  in  tens>  <five  degrees  and  hundreds  of  speed> 

<speed>  I  /////; 
<wind  direction  in  tens>  ::=  <digit>  <digit>; 
<five  degrees  and  hundreds  of  speed>  ::=  <digit>; 
<speed>  : : =  <digit>. 

7.   Metalinguistic  Symbols 

Diamond  brackets  group  the  descriptions  of  non-terminal  symbols.   The 
::=  is  the  replacement  operator  and  j  indicates  "or".   Semicolons  indicate  the 
end  of  a  metalinguistic  statement  and  periods,  of  a  section. 
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%  TEMPPOlNTtR  PLACro  DISTANCE 


ooor?oo 
0001 ?io 
ooo<  ??o 

000i?30 
0001 ?flC 
000i.?50 
000l?f0 
000?.?70 

oooi^eo 

0001P90 
0001 ?00 
0001.  ?10 
0001 3?0 
0001.  ?30 

S  0001 ?flO 
0001 350 
0001 360 
0001 370 
0001.380 
000090 
0001 flOO 
0001 fll 0 
O0OM?0 
0001.^30 
0001.  flflO 
0001.  fl50 

EDOOOi ^60 
OOOiflTO 
000i.fl80 
000ifl90 
O00i«i00 
0001510 
000l5?0 

oooi'sao 

000l5fl0 
0001550 
OOO1560 
0001.570 
0001580 
0001590 
OOOllSOO 
OOOliSlO 
OCOi.6?0 
0001630 
000^40 
0001650 
0001660 
0001670 
0001680 
0001690 
0001700 
0001710 
OOO1770 
OOO1730 
0001740 
0001750 

OOO1760 
0001770 
0001780 
0001 790 
OOOlflOO 
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t    pFTlTRMlNjr    Tr    STATION    IS    WITNTW 
I  01  r  .Rpr ^     .RMfnsoi 
jMi  fARRnCI 

fCTr  J    Aigr    n 

<F Tr 1  M    AMP    Fl > 

T     ir     STATION     IS    WTTtJlM    RMID*     PI  AC 
I  DA  ►'TPMTI 

I  01  f  ,C^iTP),ArrwiOHT| 

S'.  TTf.rlAfO-TTwppnTNTrPI 

rxr^^L{.rLA^'^SI^R» 
^TA  i-lnwTI 

1 OA  fSI 


I  nFTFRMTvr  IF  ST 
I  Dl  (  .BHf 
iMl 
I  OA 

'.01  (.TNT 
SLTTf  .n 
rKfhl  f  ,r 

»  PfTfRHTsT  PADTU 

ADM 

'.  0L(  .PPf*  > 

^^TF 
S'^TFl 

^Tfl 

^TTT 

S^^TFl 

ADM 

AOM 

LOl  r  .ODC) 

J  Mr, 

STTF 

s?:Tr 

LDI  f 

L0« 

STA 

SETF 

SfTF 

LOI  ( 

LDA 

STA 

t  INSTRT  ALL 
CLRA 
LOP 
SLIT 
rXCA 

t    1    TNSFRT    PP 
f    OFTE:PMT^F     I 

LDL( 
JA" 
^ETF 
SETFI 
SLIT( .FL 
rXCMi  (  .F 


1 
.ROC 


1 

.fine 


0    PP 

> 

(  .FL 
Lf  .F 
INTF 
F  ST 
.BOC 


ATION 

^  .RliP 

TAR 

IIPH 

p).4n 

A  fO  «  T 

I  A  ro  s 

S    OF 

I  nw 

'  PW 

•  VL 
CAR 
T     A 

r,  A 
."L 

r  ftp 
pan 

-F 

-Fl 

MTO 
^10 

,MM 

CAR 
J  A 
b>     A 

5     .RW 

f  RO 
PAP 
-I 
-w 

)  .RM 
CAR 
PAD 

TNTTR 

tA 


IS    WITMUi 
PSO> 
RDC  i 

t; 

FtiPHTl 
FMPPniNTFPI 

icr; 
iK'fLUfNcr  r 

pnyw  J 
WT» 

Ml  ) 

one  I 

vo  rii 

''Wl 

T"S; 
Avn  -n 

rpth; 

WT) 

Ml  J 

pre  ) 
k'R  F  I 
VP  Fll 

IDI 

CI 

lUSJ 

AWO  - 

A^'P  - 

AXJ 

PPCI 

lus; 

INTO 


Fl 
Fll 


%  SOliARFO  IN  RGS 

PMIO  OF  POINT 

t  RMin*?|ScAtr  i/?^6 

f  IS  STATION  WTTHIW  MlOPLF  RADlU 

t  OISAPLF  TMOSr  OiiTSIPr  OR  RMID 

F  POINTFR  IM  TFMdOrARY  APPaY 

T   fiNP  INOICrS  FOR  TFMPORARY  AR- 
t  RAv 

t    AOpPFSS  OF  irMPnpARv  ARRAY 

t    JUMP  TO  tEmPotnTPR 
%    STOPr  (iPnATFn  PniMTFe 
f    TFmPPOINtEr  PLACrn  PtSTanCF 
f  SOUARFP  IN  RC^ 

PMAX  OF  POIMT 


«  iNPlCrs  TO  TEMPORARY  ARRAY 
t     AOPRFSS  OF  TrMPORARY  ARRAY 

f  JUMP  TO  TEmPPOINTFR 

OR  THIS  POINT 

%  TOTAL  NU^'BFR  METCTH  PFPORTS 

r  IN  RLOW 

If  K'lJMPFR    FOR    RLOW    Mjvil^    OK'E 

f  TrsT  FOR  *JL  or  CRFATFR 

%  DISARLF  fOp  RAOIus  STORF  PES 

»      NOT  IN  Plow 


»  TOTAl     NllMBrp    HFICTH    PfPORTS 

t  IN    RMIP 

t  NUMPFR  for  RMID  MINIK  ONF 

t  TFST  FOR  NM  OR  nPFATFR 

t  PiSARLf  FOR  RADIUS  STORE  PES 

t         NOT  IN  RMin 


»  REMAININC,  PES  STORE  PMAy 


AG)"FOPMPOINTFRI 
L  A  CO  $  I  C  R I 
P  INOFy  TO  RC.R 
ATION  IS  WITHIN 
)  .RLOWSOI 

CARROC) 

J  AMP  If 

y     AND    EJ 
AG)»FPRMP0INTERI 
I  AOSICRI 


STAPERPNT  TO 
t 


PLOW  OF 

f 
t 

% 

f 


MARK  START  OF  HFIGHT  ONLY 

ALL  0  POINTER 

TRANSFER  To  FQRMPOINTER 

POINT 

LOWWFST  RADIUS  OF  INFLUENCE 

COMPARED  WITH  DISTANCE 
DISABLE  IF  STATION  OUTSIDE  OF 

RLOW 


T  TRANSFER 


OOO'PIO 

000< A?0 

000'ft30 

SP00^.«40 

000^  e'^o 

000'  A^.O 
'^OO*  «70 
000<  «P0 
000!i.«9C 
000<.900 
'^00 '  91  0 
000*  9?0 

opo<9:^o 

000<  940 
0001950 
00019C.O 
000<970 

ooo'.9eo 

0001990 
000?000 

ooopoio 

000?0?0 
000^030 
000?040 
000?050 
OO0?0fO 
000?C70 
000?080 
000?090 
000?100 
000?1 10 

poo?i?o 

000?' 30 
000?1 40 
0007i?0 
000?' 60 
000?' 70 
000?' ^0 

000?<90 
000??00 
000??10 
000???0 
000??30 
000??40 
000??50 
000??60 
000??70 
000??80 
000??90 
000?300 
000P310 
000?3?0 
0002330 
0002340 
0002350 
0002360 
0002370 
0002380 
0002390 
0002400 
0002410 
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%    UPDATE 


LOS 

NUMPFR  nr 

LOK  .BDC) 

LOA 

ADM 

STA 


f  A) 
PFPnPTS 

.PAIRl I 
rAPRPC) 
I OWWMDJ 

i.nwwND; 


IN  LOW 


RADIUS 


SAVT  OIST^NCF  SOHARTn 


T  DETFRMI 


»  CALL  MA 


t  nrirRMT 


X  CAl  L  MA 


%    DTTFRMT 


1!  CALL  MA 


TAFTr 
IPFTlJ 


SUP 

NFy 

FNT 


P  AL 
RK  T 


POUT 
TSMJ 
F  AT 
RY  ! 


PI  I 


MF  THF 

SFTF 

SETFl 

I  DA 

ADM 

SFTr 

SFTFl 

^FTF 

ADM 
I  OFri 
STA 

ycniiNT 
SLTK.p 
SLTTC.F 
FXCMLr'. 
MF  »'AXT 
L  OA 

SFTf 
SFTF  t 

«0M 
«;rTr 

AOM 
I  OFf  1 
STA 

XCniiNT 
SLTK.P 
SL!T( .F 
FXfMLf , 
MF  MAX! 
LDA 
SETT 

SETFl 
ADM 
SFTr 
SFTf  1 

AOM 
LOFfl 

yCO'iNT 

S'-TTf  .P 
S'-ITC.F 
FXrwLf . 
I  P'"TNT 
n  CALLF 
I  0'  (  .F[ 
S  T  L  f  .  f  I, 
TN*"  fTR 
r  T  A^o 

E^jTPY 
S  MAOr 
CSPf  .fl 

SETT 


or 


PFP0RT5  WITH  BOTH  HEIGHT  And  wINO 
%    ENARLF  ALL  PFS 


LOWBOTH 

niSAPLE 

fOr  plow  points 

oisarle 

fOR  RMin  POINTS 

MAX  FOR 


ARRAY  nr  nrboth 

t    PASS  ARRAY  ADDRESS 


COUNT 


REPORTS 


El 

Ell 

Zt 
Ell 


FOR 


TOTAL  NLiMREP 
F  OP  -El 
F  OP  -El 
fPI 
MIOPOTHI 

-I  A  Kin  El 

-r,  AND  Ell 

-J  AND  El 

- w  A  V  n  El; 

HPRflT^^I 
CARAl L) 

M  p  p  n  T  u ; 

TO  nFTFRMTNF 

or  )«MAxpnTHi 

LAG)»MAXCOUNTI 
F  l_  A  CO  S  I  r  P I 
•^ilM    ^UMRFR    Pf 
'  OWHTI 

-  I     A  M  0 

-  r,     ft  M  n 
UIDMTI 

-  J    A  K'O 
-I-'     t\"^ 
lifWTI 
r  A  R  A  t  LI 
VPHFIGHTl 

TO    f^FTFPMINF    MAX 

DC)»MAXMriGHTI 

LAG)«MAXCOUNTI 

FLAfOSTCPI 

MilM    NUMPFR    PF 

L  nwwTNDl 

-I  A  NO  El 

--"■  ANO  Ell 

t'lDwiVOI 

-J  AMO  El 

-w  AV^  Fll 

HPWT».OI 

**  A  R  A  L  L  I 
Tn  nrTFPMjsp 

nc  )«MAywlNOi 

L  AG>«MAyCOUNTI 

p"  L  A  ro  s  T  r  p ; 

S  HAVE  PFFN  CMFCKEO 
R 

AC,  ).  LINK  I 
A  r, )  »  T  C  P  I 

MPOTNTFP  RUTLOS  STAPFPPNT  BY  U5F  OF  POINTERS 
►"■XTSTOPF.  THF  MFW  POINTER  MUST  PF  IN  RGR, 
IS  NOT  CHANGED, 

PY  Fycwi  PN  .FLAG,  THFREFORE » . FL A G  PIT 
»r,).  INMTRWOI  t    MARK  START  OF 

f    INNER  WOPD 
F  OP  '£1  %    I  PIT  SET  FOR 


1!  TRANSFER  TO  MAX 
WITH  HEIGHT  ONLY 
%    LOWHT 
t    ni?APLF  FOR  PLOW  POINTS 

%    DISABLE  FOR  RMIO  POINTS 


«  FNARl F  ALL  PES 

ARRAY  OF  NRHElGHT 

1!  PASS  ARRAY  ADDRESS 


REPORTS 


%    TRANSFER  To 
WITH  WIND  ONLY 


MAX  COuNt 


%    DISABLE  EDR  PLOW  POInTS 


%    DISABLE  EOR  RMID  points 


MAX  FOP  ARRAY  OF  NRWIND 
J!  PASS  ARRAY 


t    TRANSFER  To  MAX 
FOR  ALL  STATIONS. 


ADDRESS 

COUNT 


0  MUST  BE  RESET, 
PROCESSING  ON 

INNER  WORD 
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FILI 

rpPMpniMrRi 

LOP 

LPP 
IDB 

I  DA 


P?l 


«;wAPA  ) 
i.or 

I  OP 
I  Oy 

SMAR 

LD«^ 

LDA 

LDL ( .BDC) 

AMD 
PTAl 


SETT 

PR 

STA 

I  OF 
JXI  n 
SETT 

LOX   -561 
LDF 

LOR 
LDX 
LD5 
SETF 

LOT 

LOP 

LOA 

SHAL 

OR 

STTF 

SHAPAi 

SETF 

STA 

CTSBFf  .FLA 
CRPr .FLA5) 
LDP 
LOO 

LOA 

LDl  (  .FLAG) 
t  SUPROUTINF  PuILPS 

riiL 

TEMPpOIKTrPl 

LOR        SRI 

LD8        tOJ 


»P?I 
71 

Til 
SRI 

SDI 

r  AND  If 
F  AND  -El 

NFXTSTGRFJ 
T  AND  EU 


.AI.LI 

SAI 
SPI 
161 
SAI 

.«p; 

.MASKOUTOini 

.PDCI 

♦  -?! 

G  AND  El 

tSTAPERPNT) 
tSTAPERPNTi 


.ALI  1 
J  AND 


.ALLJ 

SXJ 

SSI 

J  AND 

•  i; 

.ALLI 

SBI 

161 
SSI 

!  AND 


e; 


El 


El 


f,  AND  Zf 
NEXTSTORFI 

fO.INNERV*D»F 
. INNERWDI 

SBI 
Til 

SICRI 

rULI  WORD  P 
?l 


X  SKIP  TNlTlALTZATTON 

f   ALir.N  ON  Even  svll ari  e 

«  SAVE  DiSTANCr  IN  T1 

%  SAVE  MPOE  PATTERN  wTtM  TNOfX 

i;    IN  rgR 

%  SAVE  r  PIT  IN  G  BIT 

f  I  PIT  PESEt  FDR  noTEP  WORD 

t  DFTTRMINE  LflCATlPN  FPR  POINTER 

*  enarl*"  when  ppocrssiwr.  inner 
t   wnpr) 

»  POINTER  TO  OUTER  WOR'> 

T  BFGiN  PROCESSING  OF  INDEX  IN 

t  nilTFR  WORD 

%  SWIFT  PART  TO  RGy 

t  FIND  INDEX  PART 

r  ANP  place  in  PCS  FOR  FETCH 

r  INDEX  FOP  CURRENT  jtlTATION 

51       t  MASK  OUT  MODE  BITS 

r    SHIFT  INTO  PPINTTR  POSITION 

%  WITH  TWO  LOW  BITS  PROPPED 

t    RESTORE  E  BIT 

%    CONCATENATE 

f    STORF  POINTER  FOR  US?  U'  ANALY' 

%  SIS  ROUTINT 

%    RESUME  PROCF^SINft  OF  TNOEX 

j:  COMPUTE  nExT  SHIFT  PART 

«  FNAPLE  only  pes  CnMPiETiNG 

%         A  KORD  Fpp  STAPERPnT 

%    SHIFT  PART  TO  5(^ 

t    INTERCHANGE  RGX  S  rG^ 
X    INDEX  PART  TO  RQX 
X    SHIFT  PART  TO  RGS 

f    AND  INCREMENT  INOEX  PART 
%    IN  ORDER  To  SAVE  MODF  PATTERN 
»  INDEX  PART  TO  RGB 
»  AND  INTO  RGA 

t  INTO  PPlNTFR  WORD  POSITION 

%    POINTER  FORMED  IN  OUTER  WORD 
t    DISABLE  WHEN  PROCESSING  OUTER 
%         WORD 

r    OUTER  WORD  DUPLICATED 
%    RESTORE  F  PIT 
t    CN| Y  APPROPRIATE  HALF  WORD 
f  STORED 

ORMPOlNTFR-1 I  *  SKIP  TO  00  INNFR  WORD 
t     INNER  WORD  JLiST  COMPLETED 

t    RESTORE  MODE  PATTERN 
T  RESTORE  RGA 

r    RETURN  TO  CALLING  SEoUENCE 
OTNTERS  TO  TEMPORARY  ARRAYS. 


f  SAVE  MODE  BITS 


OC03030 
000?'^40 
0003050 
0003060 
0003070 

0003080 
0003090 
0003100 
00031 10 
00031P0 
00031 30 
00031 10 
0003150 

0003^60 
00031T0 
0003'.eO 
0003190 

ooo3?no 

0003?10 
0003??0 
0003?30 
0003?A0 
0003^50 
0003?<^0 
0003?70 
0003?80 
•0003?90 
0003300 
0003310 
0003320 
0003330 
0003340 
0003350 
00033AO 
0003370 
0003380 
0003390 
0003«00 
0003410 
0003420 
0003430 
0003440 
0003450 
0003460 
0003470 
0003480 
0003490 
0003500 
0003510 
0003520 
0003530 
0003540 
0003550 
0003560 
0003570 

0003580 
0003590 

0003600 
0003610 

0003620 
0003630 
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LDR 
LOB 

LDX 

LOfFl 

SWAPA) 

LOP 

LD*. 

LDA 

LDP 

I  on 

SETFl 
STP 

LOB 

I  00 
<;ETr 

STP 

I  orri 

I  DP 
LO? 

I  o« 
^wapa; 

LDR 
LOA 
LOP 

ion 

LCI  ( 
ftOW 
STLC 
%    STORAGF  LOC 

MAySTAPFRPt  I 
POTNTOROrpi    R[ K 

ICllRRFNTi  8LX 

JCIIRrFNTi  Rl.K 

KiFyT^TCRTi  RLK 

KrvT^MlFTi  PL»< 

K'OPPTHl  RLK 

KMjMpfRLOWj  RL*^ 

wTPPrlNTrRl  PLK 

(iPPOTKIPRi  oiv 

I^TAPrRPNTi  0L'< 

PI  fin  J  Qi_K 

Tit  RLK 

T?l  BL*^ 

T  3  !  R  L  •< 

It;TARTl  Bi_K 

jS^APTj  PL'' 

►'inSTATlnNSI  bl'' 

L'P^'TATinwSl  BLK 

^'RPOTHi  RtK 


SAI 

.ALU 

tSI 
SAf 

sp; 

*R) 

Fl  AND  -FU 

*0(.CNTR)J 

«P; 

f  PI 

F  ANO  -F; 

fOC.CNTRJJ 

.ALL) 

^Pl 


•El 


F  AND 

*Pi 

»R| 

.prC5  .PAIRl) 

.PDC) 
,Fl  AOSICRI 
ATTDNS 

FOU  2?0t 

?J 

i; 
i; 

II 
II 
'I 

i; 
?. ; 

'?8| 
flfll 

<  I 

Tl 
1  > 
1  > 

111 

•51 

•  I 


«  OUTFR  INDEX  TO  RGX 

T  FNABLF  POInTFR  FFTCH  AND  SwAPA 

*  INNFP  INDEX  TO  OUTER  WORD 

i:  DISTANCE  SQUARED  OUT  OF  RGS 

%  INNER  INDEX  TO  RCR 

»  RESTORE  Mopr  BITS 

*  disaplf  INnFR  WORO 
i;  STORE  outEr  word 


It  DISABLE  OUTER  WORD 

%  STnRF  INNER  WROD 

%  INNER  INDEX  TO  RRA 

*  DISTANCE  Squared  to  rgs 

r  innFp  index  rack  to  inner  word 

«  DISAPLE  INNER  WORD 

%  OUTER  INDEX  PACK 
%         INTO  "GA  OUTER  WORn 

f  RESTORF  MODE 


%  INCRFMFNT  FOR  STORF 

%  RETURN 

f  MAXIMUM  STAIONS  USED  IN  ANY  PE 

f  EACH  HALF  BYTE  PniNTS  TO  A  GRI 

f  POINT 

r  FOP  TEMPORARY  STORAGF 

«  POTNTFR  TO  STAPFRpNT 


%  POINTER  TO  MIDSTATONS 

y  POINTTR  TO  UPSTATIONS 

t  POINTTPS  To  STATIONS  TO  RE  USE 

t  IN  ANALYSIS  AT  FACM  POINT 

%  ARRAv  rOR  STORAGE  OF  FiFLO 

%  BEING  ANALYZED 

r  TEMPORARY  STORAGE  LOCATIONS 


V  R  M  T  I  H  I.  K   1  ) 

•^PWNDI         BLK   1| 

rilL  2*>6t 


TMr  IFF 
THF  I nw 
FOP  TFM 
DFTrf! 
NUMBER 

MEIGH 
NUMBrp 
ONLY 
I  NUMBER  0 
It    WIND 
«  5TA0ATA 
X  ROW  B 

r  DROPP 


T-MOST  COLUMN  IN  BLOCK 
ER  ROW  IM  BLOCK 
PORARY  STORAGr  WHILE 
MINING  STATIONS  USED 
OF  RFPORTS  WITH  BOTH 
T  AND  WIND 
OF  RFPORTS  OF  HEIGHT 

F  RfPoRTS  WITH  ONLY 

MUST  BE  ALIGvEO  ON  FOU 
OUNDAPY--LAST  TWO  BITS 
ED  FOR  STAPOInTER 


000  3 
0003 
0003 
0003 
0003 
0003 
0003 
0003 
0003 
0003 
0003 
0003 
0003 
0003 
0003 
0003 
0003 
0003 
0003 

0003 
0003 
0003 
000  3 
0003 
0003 
0003 
0003 
0003 
0003 
0003 

0003 
00003 

0003 
000  3 
0003 
000  3 
0004 

OOOfl 

OOOa 

OOOa 
OOOa 
DOOOfl 
OOOa 
OOOfl 
OOOfl 
OOOfl 
000-9 
OOOA 

oooa 

OOOa 
OOOfl 
OOOfl 
OOOfl 
0004 

OOOa 
oooa 

000/i 

0004 

ROOOA 

oooa 
OOOa 


(S«0 
«50 
660 
670 
680 
690 
700 
710 
720 
730 
740 
750 
760 
770 
780 
790 
800 

eio 

620 

830 
640 
6S0 
660 
870 
680 
690 
900 
910 
920 
930 
940 
950 
960 
970 
980 
990 
000 

010 
020 
030 
040 
050 
060 
070 
080 
090 
100 
1  10 
120 
130 
140 
<50 
160 
170 
180 
190 
200 
210 
220 
230 
240 
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STADATAl 


RLX   axMAXSTAPERPTI 


f  fU  SYMPpL  PEriMTlONS 

.HPCi  rou  SCOI 

.ALL  :  rou  »C1 I 

.c»jtri  rou  SC2I 

.rL*r,i  Fou  SC3I 

,STrP«U'^TIL256|  reu  SDOl 

.MASxpUTOTOPPi  ''PU  *0U 

.MAs^OL'TOTOi^i  rou  *o?; 

.STroMflUKiTlLOJ  fcu  «03| 

.PA  I*?"  I  «"0u  *r4j 

.CI nv«oi  rou  *D5^ 

.Biipe^oi  «"0U  »C*8; 

.  o  «  T  cj  5  I  «■  I?  U  S  0  9 .! 

.ftnri'Pt  f/jij  Sn^oi 

.anrvipi  rQ(_i  *Dlij 

.MixvORIlPi  «"0U  »D6?j 

.STn«»*"STrPMA  •  «'0U  *C59; 

.TrMDi  rou  *D6P; 

,f  »jTP?;TPPr  I  fQi<  »P6i; 

.MAyvppMyoi  rou  S063j 
FSf^. 


»  STPPAGF  FOR  PATA  TO  "E 
%         FACH  Pr  FOR  ANALYSIS 


f  CU  FLAG  WORD 


USFD  IN 


000A?50 
000a?<S0 
000A?70 

oooA?eo 

000/»?90 
0004300 
0004310 
OO0A3PO 

0004:^30 

000A340 
0004350 
0004360 
0004370 
0004360 
0004390 
0004400 
0004410 
00044?o 
0004430 
0004440 
0004450 
0004460 
0004470 
0004480 
0004490 


'^h 


ROUTINE  FIELDSCAN 
P  E  G  I  N 

»  SCANFIFLD  PERFOP 
»  DEfINF  statemfnt 


SCAMFlELPi  FILL 

SLIK.PDC 
LOADC .pdc 

%  INITIAL IZF  FOR  F 
IDI  (  ,CWTR 
LDA 

SLiT(.pnr 

*  FIND  PniNTFR  TO 
SFIl      SNAP 

ALIT(.PDC 

LDX 

LDA 

STA 

STL(.RnC) 
%    FIND  INDEX  TO  FT 

LOS 

LDL( .BPC) 

!ML 

SETF 

SETFl 

LDLC.PnC) 

ADM 

SETF 

SETEl 

LDL(.BDC) 

ADM 

?;ETF 

!;ETF1 

STA 

LDA 

LDEFl 
%    FIND  FIRST  GUESS 

LDB 

LDX 

SWAPAI 

LDS 

LOS 

LOA 

SETF 

LDA 

StTF 

STA 
LDB 

LDX 

SWAPAl 

LDS 
LOS 

LDA 

SETF 

LDA 

SETF 

STA 


MS  TH 
S  FOR 
DEFI 
PEFI 
nEFI 
?l 
)»MAX 
),MAX 
IPST 
).STF 
POIN 
)«?8l 
NEXT 

0(  ,B 
)»-4j 
CARC 
*RAn 
PADI 
,PTS 
ELD  F 
SAj 
.PAI 
CAR 
I  AN 
fi  AN 
.PAI 
CARP 
-E  A 
-El 
.PAI 
CARB 
E  OR 
Fl  n 
FIEL 
SRI 

CARA 

VALU 

SAI 

SBI 


E  AOJU 
CU  RE 

ne  car 

NE  CAR 
NE  CAR 

BOTHST 

cycle; 

POINT 

PlUNTI 

TOROER 

POINT 


STMENT  OF  THE  FIELD 

GISTERS 

BDC««COff I 

ALL"SC1##I 

CNTR«*C2#f; 

I!  ONTO  A  WORD  ROUNPARY 
AJ 

t    START  FETCh  OF  MAX  CYCLES/POl 


L15I 
I 


ntr; 
in 
us  J 

MIFTi 

OR  THIS  POINT 

R8; 
PDCJ 

D  e; 

D  Ell 

R28j 

nCl 

ND  '1} 

AND  -Ell 

P33j 

nci 

-El 
R  -Ell 
DINDEXI 

LLI 
ES 


SAI 

SPI 

♦FIELDI 
F  AND  -El 
fFIELOl 
I  OR  El 

GPFIELDI 

SAI 
SPI 

SAI 

SRI 

•TEMPI 

r  AND  -El 

fTEMPI 

T  OR  El 

f.PTEMPi 


t   SHIFT  TO  Unpack  first  point 


f  RADII  STORED  IN  SAME  ORDER 
»    AS  POINTS  ARE  USED 


t    TEST  FOR  ROW  OF  POInT 


t    INDfV  FOR  ROW  1  PDINTS 
r    FOR  POINTS  ON  SECOND  ROw 


%  INNER  INDEX  TO  OUTER  WORD 

%  AND  TO  RRS 

%  FETCH  OUTER  WORD 

f  AND  PROTECT 

f  FETCH  INNER  WORD 

I!  RESTORE  MODE 

»  TEMPORARY  STORAGE 


«  INNER  INDEX  TO  OUTER  WORD 

%    AND  TO  RgS 
I  FETCH  nuTFR  WORD 

%    AND  PROTECT 
r    FETCH  INNER  WORD 

%    RESTORE  MODE 
%    TEMPORARY  STORAGE 


0000010 
0000020 
0000030 
OOOOOAO 
0000050 
0000060 
0000070 
0000080 
NT0000090 
0000100 
0000110 
0000120 
0000130 
OOOOMO 
0000150 
0000160 
0000170 
0000160 
0000190 
0000200 
0000210 
0000220 
0000230 
0000240 
0000250 
0000260 
0000270 
0000280 
0000290 
0000300 
0000310 
0000320 
0000330 
0000340 
0000350 
0000360 

0000370 
0000380 
0000390 
0000400 
0000410 
0000420 
0000430 
0000440 
0000450 
0000460 
0000470 

0000480 
0000490 
0000500 
0000510 

0000520 
0000530 
0000540 
0000550 
0000560 
0000570 
0000580 
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IKTTIALIZF  TrMPPRARY 
fLR*l 
STA 


^TnRAGF  LOCATlnNS 


MEXTSTORFI 
STA         WTSUMBOTHj 
5TA        WTSllMW^Ol 
STA         WTSlJI^HT) 
STA         WTSU^'TEMP; 
STA         MRHTI 
S^A        K'RBnTW) 
STA         ^RWMOI 
STA         VRTrMP) 
»  FOR^  CnnNTFR  FPR  STATIONS  OF  TYPF 
SLTT(  ,CMTR)Bl  I 
f Shi ( .C^TR)?«| 
rArn(  .f^'TR).MAXCYCLE» 

rswi  ( .f  N'TR)?a  J 

SLIT(.CN'TR)«1) 
»     CPmPLFTF     PFM     INTTIALI7ATI0N    FOR    SCAN 

CTSPFr  .ri  ArO.SCAN4,SF?j 

STA  wTSliMPOTM) 

STA  WTSlJMWNDl 

t    thr    FTP    sca^s    ?»3»a 
«     TNTTIAl TZF    vrxTSHIFT 

I  01.  (  .RPf  >     .PAIP?«I 

LOA  CARPnC) 

STA         NFXTSMIFTI 
I  rOMPljTr  cnPRFCTTOv  FOR  THIS  STATION 


BOTH 


^F?| 


SLIT(  ,Fl  AG)«UNP4CK; 
FXfhK.Fl  ArosiCPJ 
CTSPT(.FLAro.SCANl/SF«| 


LOA 

LOLf .ROC) 

SBM 

LOS 

LDA 

SBM 

MLM 

ADM 

SHAp 

ADV 

SB^ 

ML»^ 

SHAP 

aD^- 

ADM 

STA 


HOFKlNni 
.PAIR*il?SCALF> 
CARpnC I 

OrLTAH 
FA) 

vpfwind; 

f ARPPC) 
PFLTAJ) 

Al 

hfight; 

r.PFiFLD; 

wFir.wTi 

AJ 

tst 

WTCOPBOTHI 
WTCOPBOTMJ 


t  iNfPFMFKT  COUNTFRS  FOP  DENOMINATOR 
CTSPF(.FLAG),SCANA,SFI1 j 


INCRFMFNT 
LIMIT 

STARTING  VALUE 
SKIP  UNLESS  SCAN  A 


TRANSFER  TO  UNPACK 
SKIP  FOR  SCAN  1 

5l?ISCALE  1/16 

CONVERT  DISPLACEMENT  FORM 

SIGN/MAGNITUDF  FORM 
I  COMPONENT  or  CHANGr»SCALE 

SCALE  1/256 


t  TO  SlfiN/MAGNTTUPE  EOpM 

i:  GADIENT  CHanGEiSCALE  1/?56 

%  SCALE  TO  1/16 

*  rOPRFCTON  COMPUTED 

r  WEIGHTED  CORRECTTPNISCALE  1 

%  SCALE  TO  1/16 

r  CORRECTION  CALCULATED 


f    SKIP  UNLESS  SCAN  4 


SFIli 


t  NOm 
SFAi 


LDA  WTSUMBOTHI 

ADM  weight; 

STA  WTSUMBOTHI 
LDL(.PDC)  .PAIRW 

LDA  NRBOTHj 

ADM  CARRDCJ 

STA  nrbpth; 

COMPUTE  CPPPECTION  ON  BASIS  OF  HEIGHT  ALONE 

LDA  height; 

SBM  aveheighti 


0000'590 
0000<S00 
0000610 
00006P0 
0000630 
OOOCS^O 
0000650 
0000660 
0000670 
0000660 
0000-^90 
0000700 
0000710 
00007?0 
0000730 
0000760 
0000750 
00007AO 
0000770 
0000760 
0000790 
0000600 
0000610 
00006?0 
0000630 
00006«0 
0000650 

0000860 
0000670 
0000680 
0000690 
0000900 
0000910 
00009?0 
0000930 
00009^0 
0000950 
0000960 
0000970 
0000960 
0000990 
OOOlOOO 
0001010 
00010?0 
/2560001030 
OOOlOftO 
0001050 
0001060 
0001070 
0001060 
0001090 
0001100 
OOOlllO 
0001120 
0001130 
OOOllAO 
0001150 
0001160 
0001170 
0001180 
0001190 


TO 
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SHA 

ADM 

^TA 
«  INCREMP^T 

CTS 

LDA 

aDV 

STA 
SFI?i     LDLf.BDC) 

LOA 

ADM 

STA 
»  COPRCTTHN 

CTS 

IDA 

MLM 
SHA 
ADM 
STA 

t   tncrfmfmT 

CTS 


FPR  T 

.PTc.r 


COUNT 


LDA 
ADM 
STA 

SFI3I  LDL 
LDA 
ADM 
STA 

SF5i  TXl 
SLI 
I  OA 
LDA 
LDL 
SBM 
JML 
SET 
SET 
LDL 
LDA 
SET 
SET 
SET 
SET 
IDA 
LDL 
ADM 
STA 
SET 
SET 

%  FORM  CnUNT 
SLI 

rSH 

CAD 
CSH 
SLI 
%  rOMPUTF  HE 
SFM  LDA 
SBH 
MLM 


WEIGHT) 
P       til 

WTCORHTj 
WTCnRHTj 

CnUNTERS  FOR  DENOMINATOR 
PF(.rLAG),SCANi»,sri?l 

WTSUMHT) 

WEIGHTI 

WTSUMHT) 

.PAIRl; 

NRHTI 

CARPnCI 

NRHTI 
FMPFRATURE 
LAG),TEMPANAL*SF5j 

TEMPFRATllREI 

AVETEMP; 

wriGHTl 

«l 

WTCORTEMPI 

WTCORTEMPJ 
FRS  FOR  DENOMINATOR 
LAG5.SCAN^»SFI3j 

WTSUMHT) 

WEIGHTJ 

WTSUMHTJ 

.PAIRlj 

NRTEMPj 

CARBPC; 

wRTFMPj 
NTP)#SF2J 
r)«MAXHTSTAJ 
O.MAXCYCLFJ 

NEXTSHIFTJ 

.PAIRBj 

CARBDCj 

"01 

J  AND  e; 

H  AND  El 
.PAIR24J 

CARRDCJ 

F  OR  "El 

E  OR  -EJ 

J  AND  El 

H  AND  Ell 

NEXTSTORFI 

.  P  A  I  R  1  I 

CARRDCI 

NEXTSTORFI 

r  OR  -El 

F  OP  -El 
FOR  HEIGHT  ONLY  STATIONS 
CNTR)»1I 
I ( ,CNTR)24| 
D(  .CNTRKMAXCYCLEI 
I (.rNTR)24| 
T(  .CNTR)»1| 
IGHT  DNIY  CORRECTIONS 

HEIGHTI 

AVEHFIGHTI 

WEIGHTI 


%    WflGHTfD  CORRECTIONISCALE 
J!  SCALF  TO  1/16 


1/2560001200 
0001210 


(  .BDC) 


TM(.C 
T(.BD 
P(.PP 

r  ,BPC) 


F 

El 

(  .BDC) 

F 

Fl 
F 
Fl 

(  ,BDC) 


F 

Fl 
FR 
T(, 


t    SKIP  UNLESS  SCAN  U 


t    SKIP  IF  NO  TEMP  ANALYSIS 


f  WEIGHTED  CORRFCTT 
t    SCALF  TO  1/16 


t    SKIP  UNLESS  SCAN  A 


NEXT  STATION 
START  FETCH  OF  MA 


DECREMENT  SHIFT  C 
TEST  FOR  COUNT  LF 
ENABLE  IF  SO 


FRS 


f  RESET  SHIFT  COUNT 

t  ENABLE  EOR  STORE 

%  ENABLE  IF  NEW  ROW 

t  FOR  UPDATE 

t  UPDATED  POINTER 

f  INCREMENT 

X  LIMIT 

t  STARTING  VALUE 


%    WEIGHTED  CORRECTIONISCALE 


0001220 
0001230 
000l?/>0 
0001250 
0001260 
0001270 
0001280 
0001290 
0001300 
0001310 
0001320 
0001330 
0001340 
0001350 
0001360 
ONISCAlF  1/2560001370 
000l3fl0 
0001390 
OOOHOO 
0001410 
0001420 

0001430 
0001440 
0001450 
0001460 

0001470 
0001480 
0001490 
0001500 
0001510 
CYCLES/POINT0001520 
0001530 
0001540 
0001550 
)  0001560 
0001570 
0001580 
0001590 

0001600 
0001610 

0001620 
0001630 
0001640 
0001650 
0001660 
0001670 
0001680 
0001690 
0001700 
0001710 
0001720 
0001730 
0001740 
0001750 
0001760 
0001770 
0001780 
OOO1790 
1/2560001800 


OUNT 
SS  THAN 
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BDC) 


TF 
FL 


SHAR 

ADM 

STA 
»  !NrRFMFM  COUKiTF 

f  T  S  P  F  (  .  r  I 

LDA 

ADM 

STA 
SF  T«  I     LDL  ( 

LDA 

ADH 

STA 
T  rrpRCTTCN  FPR 

CTSPTC 

I  OA 

SBH 

►•LM 

SHAR 

aDM 

STA 
t     TSrPEMFKiT  CPUMF 

CTSPF(.FL 

LDA 

ADM 

«;TA 

LOL  (  .BDO 

L  PA 

ADM 

sTA 

T  X  L  T  M  f  .  f  ►.' 

WIND  ONLY  0 

riSBK  .FLA 

SL IT( .BOC 

I nADf ,RDC 
PPINTFR  PAS 

LDA 

LDl  r  .Bon 

SBM 

JML 

SfTF 

SFTFl 

LDL( .BPC ) 

LOA 

SETF 

SETFl 

SETF 

SETFl 

LOA 

LDL (  .BPC) 

ADM 

STA 

SETF 
SETFl 
I  FORM  CPUNTER  FPP 
SLIT(,CNT 
CSHL(.CNT 
CADDC  .CM 

f sml( .r^T 

SLIT(.CNT 
SrSi      SLITC.FLA 

ExrwLr .ri 


SF  !5i 


SF7i 
t     SKIP 


I  **C\/l 


41 

WTCPRHTI 

WTCPRHTJ 
RS  FPR  DENOMINATOR 
Af.),SCAN4,SFHl 

>*tsumht; 
weight; 

WTSUMHT) 

.PAIRU 

NRHTJ 

TARBPCJ 

NRHTI 
MPFRATURE 
ArO.TFMPAKAl  ,SF?»| 

tempfrati'Rf; 

AVETFMPj 

WEIGHT  I 

AJ 

WTCPRTFMP) 

WTCPRTEMPj 
RS    FPR    DENOMINATOR 
AG),SrAN4#SFI5> 

WTSiiMHTI 

WFIG^JTI 

WTSIIMHT) 

.PAJBl I 

NRTrwp; 

r  ARRf^r; 

WBTEWP j 
TR)  ( SF6; 
N    SCAN     1 

G).SCAN1,C0MPLITFI 
)«MAywNnSTA) 

).MAyrYCLF) 
T    DElIMETFR 

NEXTSHIFTJ 

,PAIR8> 

r ARBDC; 

■  0) 

J  ANP  It 

H  AMP  If 

.PA!R2a j 

PARBDCI 

F  OR  -FJ 

F  OR  -EJ 

J  AND  EJ 

M  AND  El) 

^EXTSTOREJ 

.PAIRl ) 

CAPBPCI 

NEXTSTOREl 

F  OP  '£1 

F  OR  -El 

HEIGHT  ONLY  STATIONS 
R)«l  > 
R)?4| 

R).MAXCYCLEI 
P)?4| 

P)«n 

G)»UNPACK| 
A  G  )  » I  C  P  I 


»  SCALE  TO  1/16 


%    SKIP  UNLESS  SCAN  H 


t    SKIP  IF  NO  TFMP  ANALYSIS 


WEIGHTED  CoRRECTIONISCALI 
SCALE  TO  1/16 


X  SKIP  UNLESS  SCAN  A 


%    NEXT  STATION 


0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
1/2560 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
t    START  FETCH  PF  MAy  C YCLES/POl NTO 

0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 


DFCRFMENT  SHIFT  COuNT 
TEST  FOR  COUNT  LESS  THAN 
ENABLF  IF  SP 


RESET  SHIFT  COUNTERS 
ENABLE  FOR  STORE 

ENABLE  IF  NEW  ROW 
FOR  UPDATE 


%  UPDATED  POINTER 

X  INCREMENT 

1!  LIMIT 

f  STARTING  VALUE 

t  TRANSFER  To  UNPACK 


001810 
001 fl?0 
001830 
001 8^0 

001850 
001860 
001870 
001860 
001890 
001900 
001910 
001920 
001930 
001940 
001950 
001960 
001970 
001980 
001990 
002000 
002010 
002020 
002030 
0020^0 
002050 
002060 
C02070 
002080 
002090 
002100 
0021 10 
002120 
002130 
002M0 
002150 
002160 
002170 
002160 
002190 
002200 
002210 
002220 
002230 
0022«0 
002250 
002260 
002270 
002280 
002290 
002300 
002310 
002320 
002330 
002340 
002350 
002360 
002370 
002380 
002390 
002400 
002410 
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i.ol^.bpc^  ,pair512.scaie:) 

!^eM        fARBOC) 

ylv  nri,.TAii 


LOS 
LDA 
SBf' 
MLH 
ADK 
5HA 

SHA 
ADK* 
ADM 
5TA 

»  I^(■f<FMf^lT 
CT«: 
LDA 
A  DM 
STA 

SFT?!  IDL 
L  DA 
ADM 
STA 
TXl 

%  COMPUTF  TO 
IDA 
SHA 
10? 
CT? 
LDA 
LDI 
SHA 
ADM 
ADM 

FlNDDENOMt 

LDR 
CTS 
LDA 
ADM 
ADM 

LDL 
IMG 
StT 
SET 
L  DA 
SHA 
ADM 
SHA 
ADM 
LDS 
SET 
SET 
DENOM?l  LDA 
SHA 
LDS 
LDA 


COUMT 
PFC.F 


(  .PDC) 

TM(  .C 
TAL  C 

R 

PTf  .F 


(  .RDC 
P 


SAI 

VOF 
CAR 
PFL 
SSI 

Hi 

AVE 

r,PF 

WFI 
41 

SSI 
WTC 
WTC 
FRS  F' 
L  A  G  1  , 
WTS 
WFI 
WTS 
,PA 
KIRW 
CAR 

mTp)# 
ORPFC 

WTC 

31 

SA  j 
LAG), 

WTC 
)  .PC 

f>(  . 

f  SI 

WTC 


WIND) 

BDCI 

TAJ) 


HFIGHTI 

IFLDI 

GHTI 


OPWNDl 

DRWNDl 

DP  DENOMIMATOP 

SrAfv4,SFI7l 

U  M  W  M  D  I 

GMTI 

llMWNDl 

TRll 

NDI 

RDCI 

NDI 

srfl; 

TIDN 

DRHTI 


SCANl#FlKinnFNPM| 

OPW^Dl 

DNSI 

BDOI 

nRBPTHI 


SA 
P  F  (  .  F  L  A  G  ) 

MR 
MR 
MR 

.P 
CA 

I 

G 

WT 

?) 

WT 

i  I 

WT 

SA 

-E 

-E 

MR 

31 

SA 

MP 


f  .BDC) 

F 
Fl 

R 

R 


I 

,SCAM4,nEM0M2l 

BOTHI 

HTI 

WNDI 

AIR2j 

RBDC) 

AND  E) 

AND  Ell 

SUMHT) 

SUMWMDI 

CDPBOTHI 
I 

AMD  -Fl 
1  AND  -EIJ 

HTI 

I 

WNDI 


%  512ISCALF  1/16 

f  CONVERT  OISPLACFMEmT  FORM   TO 

t       sign/magnitudf  form 

%  I  component  of  CHANGrtSCALE 
«    SCALE  1/256 


X  TO  si(Sn/magnitude  form 

f  6ADIENT  CHanGEiSCaLE  1/256 

%  SCALE  to  1/16 

f  CORRFCTOM  COMPUTFO 

t  WEIGHTED  CORRECTIONISCALF  1/25 

»  SCALF  TO  1/16 

%  CORRECTION  CALCULATED 


»  SKIP  UNLESS  SCAM  4 


J!  MULTIPLY  By  A-1/8 


%    MULTIPLY  By  Bi«l/»M/ft'l/2 


%  for  special  scan  4  denominator 

%  multiply  by  1/4 

«        By  1/2 

I  FOR  DfN0M2 


0002420 
0002430 

0002440 
0002450 
0002460 

0002470 
0002480 
0002490 
0002500 
0002510 
0002520 
0002530 
0002540 
0002550 
60002560 
0002570 
0002580 
0002590 
0002600 
0002610 
0002620 
0002630 
0002640 
0002650 
0002660 
0002670 
0002680 
0002690 
0002700 
0002710 
0002720 
0002730 
0002740 
0002750 
0002760 
0002770 
0002780 
0C02790 
0002800 
0002810 
0002820 
0002830 
0002840 
0002850 
0002860 
0002870 
0002880 
0002890 
0002900 
0002910 
0002920 
0002930 
0002940 
0002950 
0002960 
0002970 
0002980 
0002990 
0003000 
0003010 
0003020 
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SHAP 
ADM 

ADM 

I  0< 

I  OA 

SNAP 

SETF 

SET 

SHAl 

ILMI 

SETF 

SETn 

nvpH 

SHAl 
ADM 
STA 
»  f rs  IfMRLRAT 
I  OA 
ILMI 
Sf  TFl 
5E:Tr 
OVRM 
SHAL 
ADM 
STA 
TXl  TM 
HAl  T) 
»  -XHACK  tAkFS 
-^T^JT^p  TO  S 
I'F  1ST,  3PP 
>vtATS  Of     wiM 
^TATIO^  CnPR 
W«(DIST 

IS  cnMPiiTFr 
sTpRrr  IN  PF 


0(  .BDOi 
KRBOTHj 

f  S) 

f  A> 
JP> 

/I  » 

F    OP 
F    OP 


-El 
-El 


URF 


"I     AND    El 
-  r,    A  K  n    F  1  I 
»SI 

PI 

FIELDI 
FIELDI 

►^TCORTEMpj 


f  NUMERATOR 

t    SCALE  FOR  DIVISION 


t  IN  ORDER  To  BE  SAMF  l^CAiF  AS 

t  NUMERATOR  IS  ZERO  WHFN  OENOMi 
t  NATOP  IS 

»    NUMFRATOR  AFTER  DlVlOr 

t  SCALE 


-G  AND  Ell 
-I  AND  El 
NRTEMPI 
PI 

TFMPFIFLDI 
TFMPFIELDI 
(  .CNTR)*SF1 I 

NFyTSTORE  ANP 
TADATA,  AND  USE 
,  AWD  ttTH  WORDS 
0»  HETCHT,  TEmp 
niMATFS  ARE  OPT 
ANf E*?-PAOlUS*? 
AND  RETURNED  IN 
I TAI  AND  PEl TA 


v   numerator  is  zero  whfn  ofnomi' 
»   nator  is 
*  scale 

t   next  point 

nEvtshift*  fxtracts  the  next  staperpnt 
s  this  pointer  to  fetch 
OF  the  station  report, 

FRATllREdF     .TEmPANAL     M 
AINED,     AND    THE    WEIGHTING 
V(DISTANrF*?*RAnUiS*2 

RG5.     THE     I, J    cOmpOWENTS 
J,    RFSPFCTIWFLY. 


UNPACKS    M»V 
SET),     THr    I 
FAfTORl 


COMP/' 
AND  J 


OF  DISTANCE  ARE 


IHF  POI^TFf'  FRHM  STAPFRPNT  IS  TFSTED 

fhTO  BOTH  HJG)  AND  EfEl)  HITS 
liKPAfKl   F  IM        ?l 


FOR  0  AND  THF  RESULT  1%    ANDED 
%     ALIGN  ON  WORD  BOUNDARY 


t  r  i»oP 


It  riKD 


POw  AnopfS*. 

"sFTF 

StTFl 

I  DA 

LDP 

I  OX 

SWAPAI 

LDS 

LDA 

SETF 

LDA 

SETF 

I  DR 
BYTE  WITHIN 

LDA 

I  OP 

LOx 

SWAPAI 

LCS 

SWAPAI 

LDI  (  .BDC) 

S8M 


IN  STAPfppNT 
F  OR  -El 
E  OR  'It 
NEXTSTOREI 
f  Al 
$P| 

*A| 

*STAPERPNTI 

F  AND  -El 

tSTAPERPNTI 

r  OP  'El 

SAI 

WORD 

NEXTSHIFTI 

?A| 

J? 

S£| 

.PAIR8I 

CARBDCI 


X  OUTfR  INDEX  TO  RGX 

»  INNER  INDEX  TO  RfiS 

f  FETCH  OUTER  WORn 

t  AND  DISABLE 

f  FETCH  INNER  WORD 

t  ENABLE  OUTER  WORD 

t  SAVE  STAPERNT  word 


r    OUTER  COUNT  TO  RGX 
t    INNER  COUNT  TO  RGS 

X  OCCREMPNT  SHIFT  COuNt 


0003030 
00030*0 

0003050 
0003060 
0003070 
0003080 
0003090 

0003100 
0003110 
0003120 
0003130 
0003140 
0003150 
0003160 
0003170 
0003180 
0003190 
0003200 
0003210 
0003220 
0003230 
00032'»0 
0003250 
0003260 
0003270 
0003280 
0003290 
0003300 
0003310 
0003320 
000  3  3  30 
0003360 
0003350 
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